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SUMMARY 

In  automated  cartography,  there  is  a desire  to  take 
advantage  of  the  much  faster  digitizing  and  plotting  speeds  of 
the  new  rastei  scanner/plotter  devices.  However,  sequential 
processing  of  raster  scanned  data  has  shown  an  extremely  poor 
fit  in  terms  of  both  the  speed  of  processing  and  the  develop- 
ment of  the  software. 

Previous  efforts  by  Goodyear  Aerospace  Corporation  (GAC) 
for  the  U.  S.  Army  Engineer  Topographic  Laboratories  (ETL)  have 
shown  that  the  STARAN*  Associative  Array  Processor  (AAP)  is 
ideally  suited  to  the  processing  of  raster  scanned  data.  As  a 
result,  the  AAP  provides  a much  faster  processing  rate  than  that 
obtained  using  a sequential  approach,  and  the  software  develop- 
ment time  is  significantly  reduced. 

The  purpose  of  this  present  effort  was  twofold:  first, 
to  remove  some  of  the  limitations  resident  in  the  previously 
developed  experimental  software  and  second,  to  perform  initial 
investigations  into  AAP  techniques  for  the  tagging  by  height  of 
contour  line  data.  The  result  has  been  the  development  of  a set 
of  semi -product! on  software  capable  of  producing  the  basic 
symbology  types  at  all  required  resolutions  and  easily  adaptable 
to  all  existing  STARAN  configurations.  This  software  can  accept 
data  from  either  of  the  two  existing  raster/scanner  devices 
currently  installed  at  DMA  facilities.  Also,  the  contour 
tagging  investigations  indicated  an  extremely  good  fit  between 
the  AAP  architecture  and  the  solution  to  the  contour  tagging 
problem,  and  have  resulted  in  a contract  being  awarded  to  GAC 
to  develop  the  required  software. 


*T.M.,  Goodyear  Aerospace  Corporation,  Akron,  Ohio  44315 
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As  a result  of  these  present  efforts,  GAC  recommends  the 
following : 

1.  That  all  existing  software  be  quickly  installed  in 
the  DMA/ETL  Digital  Image  Analysis  Laboratory  at 
Fort  Belvoir,  Virginia. 

2.  That  once  the  software  is  installed,  the  Input/Output 
modules  be  optimized  to  maximize  the  data  processing 
rate. 

3.  That  software  be  developed  to  automatically  generate 
the  remaining  symbols  required  to  provide  an  automated 
cartographic  production  capability. 


PREFACE 


This  final  technical  report  records  efforts  and  achieve- 
ments under  the  Associative  Array  Processing  of  Raster  Scanned 
Data  for  Automated  Cartography  program.  This  program  was  conducted 
by  Goodyear  Aerospace  Corporation  (GAC),  Akron,  Ohio  and  submitted 
by  GAC  as  GER- 16327. 

The  program  was  conducted  for  the  U.  S.  Army  Engineer 
Topographic  Laboratories,  Topographic  Developments  Laboratory, 
Mapping  Developments  Division  (ETD-TD-M)  under  contract  DAAG53- 
76-C-0146.  Mr.  R.  A.  Clark  served  as  the  Contracting  Officer's 
Representative  and  provided  valuable  assistance  in  reaching  the 
contractual  objectives.  This  effort  was  implemented  by  N.  J.  Adams, 
(project  engineer),  J.  M.  Vocar,  and  K.  Losch  with  technical 
contributions  provided  by  R.  W.  Messner,  R.  W.  Lott,  and  R.  Faiss. 


TABLE  OF  CONTENTS 


SUMMARY iii 

PREFACE V 

LIST  OF  ILLUSTRATIONS vi  i i 

Section  Title  Page 

1 INTRODUCTION  1 

1.1  Background 1 

1.2  Purpose 1 

1.3  Report  Organization  2 

2 INVESTIGATION  3 

2.1  General 3 

2.2  Approach 3 

2.3  Results 7 

2.4  Recommendations 17 

3 DISCUSSION  18 

3.1  General  .....  18 

3.2  Fortran  Preprocessing  Routine  31 

3.3  Executive 33 

3.4  Tape/Disk  Routines  41 

3.5  Tape  I/O  Routines 45 

3.6  Disk  I/O  Routines 57 

3.7  Array  Load/Unload  Routines  59 

3.8  Line  Thinning  Routine 64 

3.9  Variable  Line  Thickening  Routines  ....  69 

3.10  Double  Line  Symbol  Routines  70 

3.11  Line  Editing  Routines 70 

3.12  Broken  Line  Symbology 72 

3.13  Railroad  Symbol  Routines  76 

3.14  Area  Fill  Routines 84 

3.15  Point  Symbology  Routines  87 

3.16  Line  Separation  by  Thickness  Routines  . . 102 

3.17  Vector!  zation 106 

3.18  Contour  Tagging  Investigation  137 

Vi 


TABLE  OF  CONTENTS  (CONT'D) 


A MAGNETIC  TAPE  FORMATS  

1.  General  

2.  ETL  Scanner  RLC  Format  . . . 

3.  ETL  Plotter  Format  

4.  DMA  Scanner/Plotter  RLC  Format 

B I/O  MERGING/SEPARATING  ROUTINES  . 

C COMTAL  DEBUG  AID  

L 

L. 

L. 


i 

I . 

i 1 

« 

t 

i . 


li  vii 

/ 


LIST  OF  ILLUSTRATIONS 


Fi qure  Title  Page 

2-1  Stand-Alone  STARAN  System  Configuration  ...  4 

2-2  Point  Symbols  (Lake  Istokpoga)  10 

2-3  2-Mil  Contour  Data  Source 11 

2-4  Line  Separation  of  2-Mil  Contours 12 

2-5  4-Mil  Woodlands/Outline  Input  3:1  13 

2-6  4-Mil  Woodlands  Outline  After  Thinning  3:1  . 13 

2- 7  4-Mil  Woodlands  Outline  After  Editing  3:1  . . 13 

3- 1  Relationship  of  Image  Area  to  Tape  Record, 

Array  Load,  and  Data  Block 20 

3-2  Array  Map 21 

3-3  Effect  of  Overlap  on  Line  Symbolization  ...  23 

3-4  16-K  Data  Buffer  Layout  for  4-Mil,  2-Mil,  and 

1-Mil  Resolution 27 

3-5  Tape  - Disk  Data  Reformatting 28 

3-6  Disk-Data  Buffer  Loading  (N-Array  Load)  ...  29 

3-7  HSDB  Preprocess  Parameters 36 

3-8  Table  Modification/Move  for  Executive  Use  . . 37 

3-9  Area  Fill  Processing  Module 39 

3-10  Top  Level  Flow  of  Tape  I/O  Routines 47 

3-11  ETLIOl  Flow  Diagram 49 

3-12  ETL  R.L.C.  Input  Top-Level  Flow  Diagram  ...  50 

3-13  Reformatting  (R.L.C.  to  Binary  Plot)  Flow 

Diagram 52 

3-14  BITTAB  Format 54 

3-15  Data  Correction  During  Array  Loading  ....  61 

3-16  Load  Array  (LDBYT$)  Listing  63 

3-17  Order  Data  (0RDER$)  Listing  65 

3-18  Dump  Array  (STBYT$)  Listing  66 

3-19  Variable  Line  Thickening  Routine  Flow  Diagram  71 

3-20  VNITAB  Before  and  After  Execution  of  CALRS.  . 75 

3-21  Find  Tick  (TCKFD)  Routine  Flow  Diagram  ...  78 


. i 


j 

: j 


Vi  i i 


LIST  OF  ILLUSTRATIONS  (CONT'D) 

Figure  Title  Page 

3-22  Create  Tick  (CRE8TK)  Routine  Flow  Diagram  . . 81 

3-23  Table  Lookup  Definition  85 

3-24  1-Mil  School  Symbol  89 

3-25  2-Mil  School  Symbol  90 

3-26  4-Mil  School  Symbol  91 

3-27  Point  Symbol  Editing/Thickening  Flow  Diagram  92 

3-28  One-Mil  Resolution  Boundary  Editing  Process  93 

3-29  Input  Format  for  Point  Symbols  95 

3-30  Point  Symbol  Read  Disk  Routine 97 

3-31  Symbol  Editing  Process  for  2-Mil,  1-Mil  Data  101 

3-32  Improve  Resolution  of  1-Mil  Symbols  Flow 

Diagram 103 

3-33  Line  Separation  Routine  Flow  Diagram  ....  104 

3-34  Detailed  Flow  Diagram  of  DTCTA 108 

3-35  End  Point  of  Line Ill 

3-36  Array  Layout  After  Tagging  112 

3-37  Array  Boundary  Tagging  113 

3-38  Detailed  Flow  Diagram  of  DTCTB 115 

3-39  Top-Level  Flow  Diagram  of  General  Table  Update 

Routine 117 

3-40  Array  Boundary  Table  Layouts  118 

3-41  Flow  Diagram  of  General  Table  Update  Routine  119 

3-42  Example  of  Use  of  Vertical  Array  Boundary  Table  122 

3-43  Array  Layout  Prior  to  Updating  the  Array  Vector 

Link  Table 124 

3-44  T/F  Combinations 125 

3-45  VLNKTAB  Updated  for  T/F  Combinations  ....  126 

3-46  Resetting  Leftmost  T/F  Flag  Array  Vectors  . . 127 

3-47  Top  Level  Diagram  of  Master  Vector  List 

Routine 128 

3-48  Flow  Diagram  for  Master  Vector  Routine  . . . 129 


t 


i X 


LIST  OF  ILLUSTRATIONS  (CONT'D) 


Figure  Title  Page 

3-49  Example  of  Linked  Array  Vectors  138 

3-50  Example  of  Closed  Loop  of  Array  Vectors  . . 139 

3-51  Generalized  Data  Flow  in  CONTAGRID  ....  140 

A-1  ETL  Run-Length-Coded  (RLC)  Format  A-2 

A-2  ETL  Plotter  Format A-3 

A-3  DMA  Run-Length-Coded  (RLC)  Format  A-4 


X 


SECTION  I - INTRODUCTION 


1.1  BACKGROUND 

In  previous  contracts  performed  by  GAC,  In  the  area  of 
raster  processing,  results  showed  that  a four-array  STARAN 
AAP  could  offer  at  least  a two  order  of  magnitude  improvement 
in  processing  time  of  raster  scanned  data  when  compared  to  a 
sequential  approach.  The  work  and  results  were  documented  in 
the  following  reports. 

1.  Interim  Technical  Report  (ETL-CR- 74- 1 ) entitled 
Associative  Array  Processing  for  Topographic  Data 
Reduction  March  1974. 

2.  Final  Technical  Report  (ETL-CR-74-20)  entitled 
Associative  Array  Processing  for  Topographic  Data 
Reduction  November  1974. 

3.  Final  Technical  Report  (ETL-0046)  entitled 
Associative  Array  Processing  of  Raster  Scanned  Data 
for  Automated  Cartography  March  1976. 

These  efforts  resulted  in  an  experimental  software  package  which 
was  capable  of  producing  the  basic  types  of  symbols,  namely, 
solid  lines,  double  lines,  broken  lines,  railroads,  area  fills, 
and  point  symbols,  at  4-mi1  resolution.  Also,  the  software  was 
capable  of  performing  both  line  separation  by  thickness  and 
raster-to- vector  conversion  under  limited  conditions. 

1.2  PURPOSE 

In  general,  the  purpose  of  this  contractual  effort  was  to 
produce  a base  for  a STARAN  AAP  semi -producti on  software  system. 
This  software  would  be  capable  of  processing  all  normal  resolu- 
tions (4-mi1,  2-mi1,  and  1-mi1),  sheet  sizes,  and  line  densities 
and  producing  all  previously  generated  symbols.  The  software 
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was  to  be  directly  adaptable  to  STARAN  stand-alone  systems  and 
easily  adaptable  to  STARAN-host  systems.  Also,  the  software 
was  to  be  capable  of  processing  data  scanned  on  both  the  ETL- 
IBM  scanner/plotter  and  the  DMA  raster  scanner/plotter  (RAPS), 
and  producing  data  which  could  be  plotted  on  both  of  these. 

In  addition,  two  other  software  programs  were  to  be  expanded 
to  handle  high  line  densities  under  similar  constraints  to  those 
described  above.  These  were: 

a)  Line  separation  by  thickness,  and 

b)  raster- to-vector  (linear)  conversion. 

Finally,  an  effort  was  initiated  to  investigate  techniques  for  the 
tagging  of  raster  scanned  contour  data  with  their  associated 
elevati ons . 

1.3  REPORT  ORGANIZATION 

Section  2 of  this  report,  INVESTIGATION  - describes  the 
approach,  the  results  derived,  and  recommendations.  Section  3, 
DISCUSSION-contains  a description  of  the  overall  processing 
sequence  followed  by  a detailed  description  of  each  routine 
divided  into  its  function,  approach  and  implementation  technique 
used. 

Appendices  A,  B,  and  C contain  detailed  descriptions  of  the 
magnetic  tape  formats  that  are  used  in  this  system,  together  with 
utility  routines  used  to  either  check  the  results  or  provide  some 
additional  aid  to  the  GAC  personnel  involved  in  the  generation  of 
the  software. 


SECTION  2 - INVESTIGATION 


2.1  GENERAL 

This  section  provides  a summary  of  the  approach  used  to 
perform  the  previously  defined  goals.  The  approach  includes  a 
description  of  the  testing  system  that  was  used  and  the  tasks 
involved.  This  section  also  documents  the  results  of  this 
contract  and  the  recommendations  for  future  work  in  this  area. 

2.2  APPROACH 

The  approach  may  be  divided  into  two  parts: 

1.  The  testing  system,  and 

2.  the  tasks  involved. 

2.2.1  Testing  System 

The  GAC  STARAN  Evaluation  and  Training  Facility  (SETF)  was 
used  as  the  testing  system  for  this  work.  A block  diagram  of 
that  part  of  the  SETF  used  to  generate  the  software  is  shown  in 
Figure  2-1.  The  tape  transports  used  to  read/write  the  overlay 
data  are  37-1/2  ips,  800  bpi  tape  transports. 

The  two  RK05  cartridge  disks  are  1.45  megabits/sec  with  a 
capacity  of  1.228  x 10®  16-bit  words,  a rotation  speed  of  1500  rpm 
and  an  average  access  time  of  50  milliseconds.  The  first  disk 
drive  DKO  contained  all  the  system  software  and  cartography  soft- 
ware in  the  form  of  linked  modules.  The  second  disk  drive  OKI 
was  reserved  for  the  reformatted  Input  (unprocessed)  and  Output 
(processed)  data,  together  with  off-loaded  tables,  and  semi- 
processed  data  generated  during  processing. 

The  Comtal  8000  series  display  is  used  as  a debug  aid  to 
display  four  adjacent  array  loads  of  data  after  any  stage  of 
processing. 


■CO 


Figure  2-1.  Stand-Alone  STARAN  System  Configuration 


2.2.2  Tasks 

In  order  to  accomplish  the  overall  software  capabilities 
described  in  Section  I,  Goodyear  Aerospace  performed  the  follow- 
ing tasks. 

0 Wrote  a Fortran  preprocessing  routine  which  generates 
all  parameters  required  by  the  executive  (controlling) 
routine;  these  parameters  are  generated  from  an  initial 
card  deck  input  by  the  user  containing  resolution, 
sheet  size,  symbol  type,  etc. 

0 Modified  the  executive  to  allow  the  generation  of  all 
currently  required  symbols  while  retaining  enough 
flexibility  to  add  new  symbol  types  as  required  in  the 
future. 

0 Developed  magnetic  tape  read/write  software  to  allow 
the  processing  of  raster  scanned  data  from  both  the 
ETL  raster  scanner/plotter  and  the  DMA  scanner/plotter. 

0 Expanded  the  line  symbol  generation,  and  line  editing 
capabilities  to  handle  1-mil,  2-mil,  and  4-mil  data. 

0 Generated  a new  area  fill  technique  which  does  not 
require  the  production  of  an  area  fill  symbol  library. 
The  resulting  routines  will  handle  all  required 
resol uti ons . 

0 Developed  a point  symbol  library  and  new  point  symbol 
generation  routines  to  minimize  I/O  and  yet  satisfy 
cartographic  standards. 

0 Expanded  the  broken  line  symbol  generation  software  to 
operate  at  all  three  resolutions. 

0 Generated  a new  railroad  symbology  routine. 

0 Expanded  the  line  separation  routines  to  operate  with 
user  selected  line  widths,  variable  within  pre- 
determined limits. 


0 Expanded  the  vectori zati on  software  to  operate  on  dense 
line  data  at  1-mil.  2-mil.  and  4-mil  resolutions. 


0 Developed  a Comtal  debug  routine  used  to  Inspect  a 
2x2  set  of  arrays  at  any  stage  of  the  raster 
processing  cycle. 

0 Performed  some  contour  tagging  investigations  which 
provided  part  of  the  basic  information  for  writing 
a proposal  which  was  accepted  called  Contour  Tagging 
and  Gridding  using  STARAN/CDC  software. 

The  approach  to  testing  the  above  software  was  to  use  digitized 
line  data  generated  from  various  source  materials  scanned  on  both 
the  ETL-IBM  raster  scanner/plotter  and  the  DMA  raster  scanner/ 
plotter  (RAPS). 

In  particular; 

0 ETL  used  their  ETL-IBM  cartographic  scanner/plotter 
to  perform  further  raster  scanning  at  Fort  Belvoir 
of  2-mil  and  1-mil  Lake  Istokpoga  source  data  (binary 
plot  format).  This  data  was  minimized  to  that  data 
that  would  be  required  to  fully  check  out  the  expanded 
STARAN  AAP  software  and  satisfy  GAC  and  ETL  personnel 
to  the  correctness  and  resolution  quality  of  the 
symbolized  results. 

0 ETL  also  generated  2-mil  and  1-mil  resolution  run- 
length  coded  data  for  testing  the  new  Input  routines. 

0 DMA  generated  2-mil,  and  1-mil  resolution  run-length- 
coded  data  from  their  raster  scanner/plotter  suitable 
for  testing  the  following: 

1)  The  new  I/O  routines,  and 

2)  the  line  separation  by  thickness. 
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0 GAC  demonstrated  the  line  symbolization,  point 

symbology,  and  area  fill  symbology  to  ETL  personnel 
during  the  contractual  effort.  While  most  of  the 
raster  processing  software  was  developed  for  the 
STARAN  AAP,  some  utility  software  was  developed  for 
the  Sigma  9.  These  utility  routines  are  documented 
In  Appendix  B of  this  report. 

The  following  AAP  processed  Information  was  sent  to  ETL 
for  plotting  on  their  scanner/plotter. 

0 Three  tapes  of  all  point  symbols  on  the  Lake  Istokpoga 
map  sheet  for  each  of  the  three  possible  resolutions. 

0 Three  line  separation  tapes  of  a 10"  x 10"  graphic 

which  contains  three  different  line  thicknesses  scanned 
at  l-m11  resolution  on  the  DMA  raster  scanner/plotter. 

The  three  tapes  correspond  to  each  of  the  three  line 
widths,  namely  4-m11s,  7-m11s,  and  16-m11s. 

0 Two  line  separation  tapes  of  a 10"  x 10"  contour 
graphic  which  contain  Index  contour  and  non-index 
contour  Information.  Again,  this  data  was  generated 
on  the  DMA  raster  scanner/plotter  at  2-m11  resolution. 

2,3  RESULTS 

In  general  terms,  the  results  of  this  contractual  work  are 
as  follows. 

1.  There  exists  a semi -producti on  software  package 
capable  of  performing  many  of  the  raster  processing 
functions  required  In  automated  cartography. 

2.  This  software  Is  capable  of  operating  on  a stand-alone 
STARAN  system. 

3.  This  software  Is  easily  adaptable  to  a host  computer 
connected  to  a one  array,  two  array,  or  four-array  STARAN. 
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4.  The  software  is  capable  of  handling  Input  data  from 
both  the  ETL-IBM  raster  scanner/plotter  and  the  DMA 
raster  scanner/plotter.* 

5.  The  software  is  capable  of  generating  Output  (processed) 
data  for  both  of  these  scanners/plotters.* 

6.  The  software  is  capable  of  processing  data  digitized 
at  all  the  normal  resolutions  (4-mil,  2-mil,  or  1-mil) 
for  standard  map  sheet  sizes  of  19"  x 22"  data  areas. 

7.  This  software  is  capable  of  processing  larger  sheet 
sizes  with  minimal  modifications. 

8.  The  software  is  capable  of  generating  all  the  basic 
symbol  types. 

9.  The  generation  of  the  different  basic  symbol  types 
have  been  demonstrated  by: 

a)  the  plotting  of  STARAN  processed  data  on  the 
ETL-IBM  scanner/plotter, 

b)  the  display  of  STARAN  processed  data  on  the 
Comtal  graphics  screen  existing  in  the  testing 
system,  and 

c)  the  hard  copy  output  of  polaroid  pictures  of 
sub-sections  of  processed  sheets  using  GAC's 
DICOMED  Recorder. 

10.  The  fast  data  processing  rate  obtained  with  the 
experimental  software  also  exists  at  the  higher 
resolutions  (2-mil,  1-mil). 

11.  The  software  can  perform  line  separation  by  thickness 
of  lines  within  certain  pre-defined  limits  and  selec- 
table at  run-time  by  the  user. 

12.  A capability  exists  to  perform  basic  automatic  line 
editing  during  processing. 

It  should  be  noted  that  crossed  input/output  formats 
(i.e.,  IBM  input,  DMA  output)  will  have  some  distortion 
due  to  1-mil  vs  25vi  SPOT  SIZE  (or  2 mil/50v,4  mil/lOOu). 
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13.  A technique  for  the  tagging  of  contour  lines  with  | 

their  elevations  has  been  devised  which  should  I 

provide  a base  for  a cost-effective  contour  tagging  | 

system  using  the  STARAN  and  a sequential  (host)  j 

computer.  I 

14.  The  current  software  has  built-in  debug  aids  allowing  I 

the  easy  maintenance  and  upgrading  of  the  software  I 

in  the  future.  These  aids  include  the  capability  | 

to  pass  data  currently  being  processed  to  a Comtal  I 

graphics  display  for  viewing  at  different  stages  of  I 

the  symbolization  cycle.  I 

In  addition  to  these  results,  several  raster  processing  ! 

functions  were  run  on  the  STARAN  and  the  results  plotted  using  a I 

DICOMED  recorder  to  provide  hard  copy  documentation.  Examples  of 
these  processes  are  given  below. 

A subsection  of  the  Lake  Istokpoga  point  symbols  are  shown 
in  Figure  2-2.  These  symbols  have  been  produced  on  the  DICOMED 
plotter  at  twice  map  size.  Notice  that  the  deformities  in  the 
rotated  symbols  such  as  churches,  schools  are  greatly  reduced  at 
the  higher  resolutions.  In  fact,  at  map  size  these  symbols  seem 
to  be  quite  satisfactory. 

A 4"  X 4"  subsection  of  the  original  10"  x 10"  2-mil  contours 
source  is  shown  in  Figure  2-3. 

Figure  2-4  shows  the  separated  index  (10  mil)  and,  non-index 
contours  (4-mil).  These  separations  constitute  part  of  the 
contour  tagging  procedure. 

Figures  2-5,  2-6,  2-7  show  a 2"  x 2"  section  of  the  woodlands 
overlay  during  processing.  These  pictures  were  produced  from 
photographs  of  the  Comtal  screen  and  are  approximately  3:1  in  ; 

magnitude.  The  first  picture  shows  the  input  data,  the  second 

I j 


9 


a.  NON-INDEX  CONTOURS 


Figure  2-4.  Line  Separation  of  2-Mil  Contours 
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Figure  2 


-5.  4-Mil  Woodlands/  Figure  2-6.  4-Mil  Woodlands 

Outline  Input  3:1  Outline  After  Thinning  3:1 


Fi gure  2- 7. 


4-Mil  Woodlands  Outline  After  Editing  3:1 


picture  shows  the  thinned  data,  and  the  third  picture  shows  the 
edited  (clipped)  data.  This  process  works  on  both  closed  and 
open-ended  line  segments. 

Table  2-1  shows  the  measured  performance  to  generate  a 
typical  symbol  on  the  Lake  Istokpoga  map  sheet  using  the  carto- 
graphic package.  The  time  to  generate  any  symbol  (area  fill, 
line  symbol  using  vectorizing,  line  symbol  without  vectorizing, 
line  separation)  is  approximately  the  same.  However,  the  time 
to  perform  the  array  processing  routines  such  as  line  thinning, 
line  editing,  etc.  do  vary  with  the  line  density  of  the  map  sheet. 

Table  2-II  shows  the  total  times,  including  I/O,  for  symbol 
generation  using  the  existing  STARAN  S-500  stand-alone  system  at 
Goodyear.  Since  the  I/O  depends  greatly  upon  the  overall  system 
configuration  and  the  speed  of  the  peripheral  devices  (disk, 
tapes,  etc.),  this  area  of  the  software  was  not  optimized.  However, 
in  a semi -producti on/producti on  system  it  is  highly  desirable  to 
minimize  this  I/O  by  using  fast  peripherals,  overlapping  I/O  and 
reformatting  where  possible.  Also,  the  use  of  double  buffering 
will  help  minimize  the  I/O  portion  of  the  total  processing  time. 

Conclusions  that  can  be  drawn  from  these  results  are  that: 

0 Particularly  at  the  higher  resolutions  it  is  necessary 
to  produce  clean,  "thinned"  line  data  prior  to  the 
generation  of  area  fill,  and  line  symbology.  The  auto- 
matic open-ended  line  clipping  routines  go  a long  way 
towards  providing  such  clean  data.  Future  efforts  in  the 
contour  tagging  area  should  further  enhance  the  area  of 
automatic  line  joining. 

0 The  resolutions  of  2-mil,  and  1-mil  seem  satisfactory  for 
the  generation  of  all  symbology. 

0 The  tape  I/O  and  disk  I/O  are  significant  parts  of  a 
semi -producti on/producti on  system  and  must  be  addressed 
accordi ngly. 
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TABLE  2-1  - AAP  PERFORMANCE  FOR  LAKE  ISTOKPOGA  MAP  SHEET  USING 
STARAN  S-500  AT  4-MIL,  2-MIL,  AND  1-MIL  RESOLUTIONS 


i 


ARRAY  I/O 

ARRAY  PROCESSING 

RESOLUTION 

(SECONDS) 

(MINUTES) 

4 

4.5 

0.42 

2 

18.0 

1.30 

1 


72.0 


4.23 


) 
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2.4  RECOMMENDATIONS 


Goodyear  Aerospace  recommends  that: 

0 All  software  be  quickly  installed  at  the  ETL  facility. 

0 Efforts  be  applied  to  developing  efficient  (fast) 
Input/Output  of  the  data. 

0 A technique  be  developed  to  perform  an  on-line  merging 
of  such  data  as  the  digitized  area  fill  symbols,  and  the 
merged  information  of  each  colored  overlay. 

0 A software  routine  be  developed  to  generate  automatically 
the  remaining  point  symbol  library  information  required 
to  produce  the  remaining  symbology. 

0 An  effort  be  made  to  continually  update  the  existing 
software  routines  with  new  techniques  that  may  be 
developed  through  investigations  such  as  contour  tagging 
and  enhanced  editing  techniques,  etc. 


SECTION  3 - DISCUSSION 


3.1  GENERAL 

3.1.1  Section  Description 

The  details  of  the  STARAN-AAP  raster  processing  software 
are  discussed  in  this  Section.  The  specific  routines  developed 
and/or  modified  under  this  contract  are  shown  in  Table  3-1  with 
their  instruction  counts. 

The  discussion  of  each  routine  is  broken  down  into  three 
subsections:  1)  FUNCTION,  2)  APPROACH,  and  3)  IMPLEMENTATION. 

The  first  subsection  describes  the  routine's  specific  function. 
The  second  subsection  includes  a description  of  the  general 
approach  taken  to  perform  the  function.  The  last  subsection 
details  the  implementation  of  that  approach.  In  some  cases  the 
reader  is  referenced  to  previous  technical  reports  which  contain 
material  relevant  to  the  discussions  herein  and  which  may  clarify 
the  routine  under  discussion. 

Table  3-1.  STARAN-AAP  RASTER  PROCESSING  SOFTWARE 


Routi ne 

mOggll 

Routi ne 

Instructions 
(No.  ) 

Executi ve 

300 

Line  Editing 

300 

Preprocess 

225 

Broken  Line  Symbols 

1100 

Tape/Disk 

275 

Railroad  Symbols 

820 

Tape  I/O 

550 

Area  Fill  Symbols 

135 

Disk  I/O 

40 

Point  Symbols 

1300 

Array  I/O 

235 

Line  Separation 

225 

Thin 

225 

Vectori zati on 

1200 

Variable  Line  Symbols 

175 

Mi  seel  1 aneous 

Comtal  Debug  Aid 

400 

HSDB  (Definition) 

1 

150 

TOTAL 

7655 

I 

J 


I 

J 


I 

£ 


1 

1 

1 
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3.1.2  Processing  Sequence 

It  is  very  important  to  understand  the  data  processing 
sequence  within  the  STARAN  AAP.  Since  data  movement  (read/write 
tape,  etc.)  is  a significant  part  of  the  overall  processing  time, 
it  becomes  clear  that  efforts  to  minimize  the  amount  of  data 
during  data  transfers  will  provide  a desirable  overall  performance 
advantage.  Also,  the  intermediate  buffer  storage  requirements 
become  unreasonably  large  when  processing  high  resolution  (2-  and 
1-mil)  date.  Hence,  the  use  of  run-length-coded  data  for 
Input/Output  and  vector! zati on  as  output  becomes  a very  desirable 
processing  capability. 

Figure  3-1  helps  to  define,  with  respect  to  the  image  area, 
a tape  record,  an  array  load,  and  a data  block.  The  ETL  scanner/ 
plotter  plot  tape  format  is  used  (in  this  example)  for  the  source 
tape  as  well  as  the  plot  tape.  The  tape  record,  as  shown  in 
Figure  3-lA,  is  eight  scan  lines  wide  and  as  long  as  one  image. 

In  other  words,  the  first  byte  of  the  tape  record  will  contain  the 
first  cell  of  8 scan  lines. 

An  array  load  refers  to  the  amount  of  data  processed  by  one 
associative  array  and  is  224  cells  wide  and  248  cells  long  (as 
shown  in  Figure  3-lB).  To  load  one  array,  it  is  necessary  to 
read  224  scan  lines  or  28  tape  records,  and  this  amount  of  data 
is  defined  as  a data  block  as  shown  in  Figure  3-lC. 

The  processing  sequence  is  as  follows: 

1.  Read  a data  block  from  the  source  tape, 

2.  Load  the  arrays,  process  the  data,  unload  the  arrays, 

3.  Repeat  step  2 until  the  entire  data  block  has  been 
processed , 

4.  Write  the  data  block  to  the  output  tape,  and, 

5.  Repeat  steps  1 through  4 for  all  data  blocks  of  the  image. 

The  processing  sequence  is  described  in  more  detail  in  Section  2.3, 
Testing  system. 
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(A)TAPE  RECORD 


(B)ARRAY  LOAD 


(C)DATA  BLDCK 


. 8 SCAN  LINES 


224  SCAN  LINES 


224  SCAN  LINES 


Figure  3-1.  Relationship  of  Image  Area  to  Tape  Record, 
Array  Load,  and  Data  Block 


3.1.3 


A diagram  of  the  array  map  is  shown  in  Figure  3-2.  Th 
solid  outside  line  is  the  array  size.  The  solid  inside  line 
shows  the  amount  of  source  data  loaded  into  the  arrays  at  one 
time.  The  dashed  line  shows  the  amount  of  data  unloaded  from 
the  arrays  after  processing  (in  this  case  32  cell  overlap). 


There  is  a variable  size  overlap  (N  cells)  between  array 
loads  and  data  blocks  to  eliminate  boundary  problems.  Each  array 
discards  the  outer  N/2  cells  as  indicated  in  Figure  3-2.  The  32- 
cell  field  on  the  left  of  the  array  is  reserved  as  working  space. 

The  4 words  at  the  top  and  bottom  of  the  array  are  working  space 
for  the  line  thinning  algorithm.  The  significance  of  the  overlap 
follows.  (However,  for  an  explanation  of  why  the  overlap  is  required 
during  processing  the  reader  is  referenced  to  Section  2.3, 
Implementation,  of  the  Final  Technical  Report  dated  November  1974 
and  titled  "Associative  Array  Processing  for  Topographic  Data 
Reduction  ( ETL-CR-74-20) . " ) 

Formerly,  the  size  of  the  overlap  was  determined  to  be  at 
least  as  large  as  the  width  of  the  widest  line  to  be  generated 
(32  cells),  or  originally  input  prior  to  thinning  and  symbolization. 
Expressed  mathematically: 

Overlap  Width  = Maximum  Line  Width. 

However,  Figure  3-3  shows  what  happens  to  a line  of  approximately  30 
mils  at  1-mil  resolution  that  is  thinned  and  rethickened  to  generate 
a Line  Symbol  of  approximately  the  same  width.  Because  all  of  the 
original  line  is  not  contained  in  block  N + 1 when  being  processed, 
the  thinned  center  line  becomes  ri ght- shi f ted  compared  with  the 
actual  center  line  obtained  in  block  N.  In  this  example,  the  result 
is  a duplication  of  data  after  symbolization.  Thus,  the  appearance 
of  a second  line  (Figure  3-3,  sheet  2).  The  solution  was  to  increase 
the  overlap  between  arrays  in  consecutive  data  blocks.  Expressed 
mathematically,  the  new  formula  is: 

Max.  Line  Width  Symbol  Generated 
Overlap  Width  = Greater  of:  or 

Max.  Width  of  Input  Line  Data. 

Based  on  the  possible  line  symbols  that  could  appear  at  1-mil 
resolution,  this  width  becomes  72  cells  for  a symbol  such  as  Dual 
Highway  #820  (see  Symbols  for  Medium  Scale  Maps). 


1 


For  this  reason,  the  software  was  modified  to  allow  the 
overlap  width  to  be  a variable  and  to  be  selected  prior  to 
processing.  Therefore,  the  user  selects  the  overlap  width 
depending  upon  resolution,  width  of  symbol  to  be  generated,  and 
maximum  thickness  of  input  line  data.  This  parameter  is  input 
as  a parameter  in  the  Fortran  Preprocessing  Routine.  This 
added  flexibility  (variable)  will  minimize  data  processing.  For 
example,  when  processing  4-mil  and  2-mil  data  the  overlap  can 
be  reduced  by  a factor  of  2 and  4. 

3.1.4  Disk  Implementation 

The  last  section  gave  a brief  description  of  the  general 
processing  sequence. 

Assuming  a large  enough  data  buffer  area  the  sequence 
would  be  as  follows: 

1.  Read  a data  block  from  the  input  tape  to  the  data 
buffer, 

2.  Load  the  arrays  from  the  data  buffer, 

3.  Process  the  arrays, 

4.  Unload  the  processed  data  from  the  arrays  to  the  data 
buffer, 

5.  Repeat  steps  2 through  4 until  the  entire  data  block 
is  processed, 

6.  Write  the  processed  data  block  to  the  output  tape, 
and 

7.  Repeat  steps  1 through  6 until  the  entire  image  is  ' 

processed. 
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However,  in  most  cases,  the  amount  of  data  in  the  data  block 
exceeds  the  size  of  the  data  buffer.  Following  is  a description 
of  the  basic  disk  implementation  technique  when  the  data  block 
(in  binary  coded  plot  format)  is  larger  than  the  data  buffers. 


1.  Read  as  many  data  block  records  from  the  input  tape 
that  can  be  contained  within  the  data  buffer.  (See 
Figure  3-4.) 

2.  Divide  each  tape  record  into  N-array  load  parts,  and 
create  a disk  record  from  the  same  section  of  all  the 
tape  records,  as  shown  in  Figure  3-5.  Store  each  of 
these  records  on  the  disk  leaving  enough  space  for 
similar  records  from  the  other  tape  records. 

3.  Repeat  steps  1 and  2 for  the  remaining  tape  records, 
filling  in  the  spaces  between  each  N-array  disk  load. 

4.  When  all  tape  records  have  been  reformatted,  load  one 
N-array  load  into  the  data  buffer  area  (see  Figure  3-6). 

5.  Load  the  arrays  from  the  data  buffer. 

6.  Process  the  arrays. 

7.  Unload  the  processed  data  from  the  arrays  to  the  data 
buffer  overwriting  the  input  data. 

8.  Repeat  steps  5 through  7 until  all  N-arrays  are  processed. 

9.  Unload  the  N-array  load  onto  the  disk. 

10.  Repeat  steps  4 through  9 until  all  N-array  loads  are 
processed. 

11.  Generate  output  records  in  data  buffer  from  disk  data. 

12.  Output  tape  records. 

13.  Repeat  steps  11  and  12  until  all  output  records/data 
blocks  written  onto  tape. 

14.  Repeat  steps  1 through  14  for  all  remaining  data  blocks 
of  the  image. 
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Figure  3-4.  16-K  Data  Buffer  Layout  for  4-Mil 

2-Mil.  and  1-Mil  Resolution 


The  number  of  4K  areas  used  is  determined  prior  to  processing 
and  depends  upon  the  amount  of  space  required  during  processing 
for  intermediate  storage  for  tables,  etc.  Also,  the  number  N 
of  arrays  to  be  loaded  in  each  iteration  also  is  constrained  in 
a similar  manner.  Both  of  these  variables  are  inputs  to  the 
preprocessing  routine. 

The  effect  of  the  disk  operation  on  the  overall  processing 
time  is  discussed  in  subsection  2.3  as  part  of  the  performance. 

3.1.5  Variable  Length  Input  Records 

For  binary  coded  data,  the  appropriate  disk  record  size 
is  known  when  the  sheet  size  and  resolution  are  given.  However, 
for  run-length-coded  data,  the  data  block  size  would  be  a function 
of  line  density.  It  was  not  considered  possible  to  segment  the 
number  of  input  records  read  into  the  data  buffer  during  processing 
into  fix  numbers.  This  was  due  to  the  varying  line  density.  The 
following  approach  was  taken  to  utilize  the  disk  I/O  software. 

When  processing  input  data  consisting  of  variable  length 
records,  as  in  the  case  of: 

0 ETL  run- length- coded  data,  and 

0 DMA  run- length-coded  data. 

The  input  records  were  transformed  into  their  scan  plot  equivalents. 
Therefore,  eight  run-length-coded  records  generated  one  binary  scan 
plot  record.  This  technique  was  considered  acceptable  because  it: 

0 Minimized  additional  software  modules,  and 

0 Could  be  used  to  overlap  certain  procedures  such  as 
reading  the  variable  length  records  while  reformatting 
others  to  the  plot  format. 


3.2  FORTRAN  PREPROCESSING  ROUTINE 


3.2.1  Function 

This  routine  reads  input  parameters  from  cards  and  generates 
a collection  of  data  items  required  by  the  Executive  to  control 
the  raster  processing  software. 

3.2.2  Approach 

The  approach  was  to  generate  a Fortran  routine  (PRPRC) 
which  was  executed  prior  to  processing  and  which  would  expand  as 
the  procedures  themselves  expanded. 

The  preprocessing  routine  performs  the  following  steps: 

1.  Reads  the  input  card  parameters. 

2.  Computes  the  array  parameters  for  subsequent  array 
processing. 

3.  Computes  number  of  data  blocks/sheet. 

4.  Computes  the  number  of  iterations/data  block. 

5.  Computes  the  number  of  32-bit  words  in  a plot  tape 
record  format. 

6.  Compute  all  the  disk  I/O  parameters. 

7.  Computes  the  number  of  cells  for  line  symbol  generation 
processing. 

8.  Identifies  what  process  is  to  be  performed  and  generates 
the  keyboard  requests  for  the  operator  to  load  the 
specific  I/O  and  Processing  modules  from  the  system 
disk  OKO. 


The  card  Input  parameters  are  shown  below. 


Parameter 


Integer  Value 
Right  Justified 


Resoluti on 


Number  of  Arrays  Used 
Number  Tape  Records/Sheet 
Number  of  Bytes/Plot  Record 
Symbol  Type* 

Line  Symbol  Thickness  to  be  Generated 
Number  of  Plot  Records  in  Array  Overlap 
Number  of  4096  Word  Buffers  Used 
Input/Output  Type** 

Line  Editing  Thickness 


*Symbol  Type  to  be  generated 


Value 

0-19 


Type 

Line  Symbol  Using  Vectorizing 
Intermittent  Streams 


Second  Class  Roads 


Trai Is 


20-39 


Area  Fill 


40-59 


60-79 


80-99 


100-119 


Line  Symbol  (Solid) 
Point  Symbols 
Railroad  Ticks 
Line  Separation 


**lnput/0utput  Tape  Format  Type 


ETL  Plot/Plot  Format 
ETL  RLC/Plot  Format 
DMA  RLC/RLC  Format 

DMA  RLC  Input/ETL  Plot  Output  Format 


In  addition  to  these  10  cards.  Line  Separation  processing 
requires  two  additional  parameters  (see  Line  Separation)  which 
define: 

a)  A flag  (1  or  0)  to  indicate  if  the  fattest  line  Is 
required  (1  = yes,  0 = no),  and 

b)  The  line  weight  desired. 

3.3  EXECUTIVE 
3.3.1  Function 

The  executive  routine  is  the  controlling  program,  all  Raster 
Processing  routines  including  I/O  routines  are  Initiated  by  this 
program. 


3.3.2 


iroach 


The  approach  was  to  utilize  the  existing  routine  written  for 
processing  4-m11  data.  (See  Report  ETL-0046  Associative  Array 
Processing  of  Raster  Scanned  Data  for  Automated  Cartography.)  This 
routine  was  expanded  to  allow  the  processing  of  varying  resolution 
data. 

3.3.3  Implementation 

The  executive  routine  (EXEC)  Is  written  in  APPLE  and  executes 
from  Bulk  Core  Memory. 


The  Executive  routine  has  been  designed  around  a hierarchical 
structure  whereby, all  routines  are  placed  on  one  of  several  level 
categories.  These  levels  are  directly  associated  with  the  general 
processing  sequence  discussed  earlier  In  this  section. 

A brief  description  of  each  of  the  four  levels  is  given 

below. 
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Level 

Function 

Examples  of  routines 

0 

These  routines  are  executed 
once  prior  to/after  per- 
forming an  overall  proce- 
dure on  an  overlay  sheet 

CELLSET:  sets  Line  Symbol 

Generation  parameters 
PREPRO:  Performs  Point  Symbol 
pre-processing 

MVECT:  Generates  master-vectors 

for  subsequent  Broken 

Line  Symbology 

1 

These  routines  are 
executed  for  each  data 
block  of  information 

TPRDF:  Read  a data  block  of 

Tape  I/O  data 

SETPTBT : Set  the  point  bits  in  a 
data  block  of  vectors. 

BLKPRO:  Do  point  symbol  process- 
ing for  a data  block 
of  points 

2 

These  routines  input/ 
output  reformatted 
data  between  the  disk 
and  a Bulk  Core  data 
buffer 

D2BI:  Send  I/P  data  to  buffer 

for  subsequent  processing 
B2D0:  0/P  processed  data  to 

disk  for  subsequent  outpu 

3 

These  routines  perform 
array  I/O  and  all 
processing  of  data 
within  the  arrays  on 
an  iteration  basis 

LDAl:  Load  arrays  with  I/P  data 

THINl:  Perform  Line  Thinning  on 

array  data 

AFILl:  Perform  Area  Fill. 

All  level  3 routines  are  executed  from  page  memory  (page  1 
or  page  2). 

The  main  steps  involved  in  the  Executive  routine  are  given 
below: 

1.  Read  the  input  parameters  generated  by  the  Fortran 
preprocessing  routine  PRPRC. 

2.  Modify  and  save  the  parameters  in  the  High  Speed  Data 
Buffer  (HSDB)  areas. 

3.  Select  all  Level  0-3  routines  required  for  this  pass. 

4.  Perform  all  data  processing  for  this  pass  for  all  data 
on  the  overlay  sheet. 


5. 


Repeat  steps  3.  and  4.  until  a11  passes  have  been 
performed . 

6.  Terminate  processing. 

These  steps  are  now  described  In  more  detail. 

The  Input  parameters  generated  by  PRPRC  are  read  from  the  disk 
file  F0R002.DAT  generated  by  PRPRC  and  located  on  the  system  disk 
DKO.ky  calling  and  executing  the  subroutine  RDISK. 

This  subroutine  reads  a multi-record  file  Into  the  High  Speed 
Data  Buffer  starting  at  location  bOB^g* 

Figure  3-7  shows  the  location  of  the  parameters  In  HSDB  after 
executing  RDISK. 

The  variables  BFADRFAD,  TRECAD(I),  1=1,8  and  ARECLOC(J) 

J=l,28  are  all  computed  In  the  PDP-11  using  Fortran  single  precision 
Integer  arithmetic.  For  example,  where  the  buffer  address  should 
be  BOOO^g  and  the  Reformat  area  base  address  should  be  FOOO^g,  they 
are  3000jg  and  7000  respectively.  Therefore,  the  Executive 
unconditionally  sets  bit  locations  0 and  16  of  all  the  above  32-b1t 
values . 

The  executive  then  stores  the  16-b1t  values  TRECAD(I)  1=1,8 
In  a 32-b1t  word  table  RECADR(I),  1=1,8.  The  executive  also  stores 
the  16-b1t  values  ARECLDC(I)  1=1,28  In  both  16-b1t  halfwords  of  the 
array  record  table  RESETLOC(I)  1=1,28.  Finally,  the  Base  Disk 
Block  numbers  DADDR(I)  1=1,36  are  saved  as  32-b1t  values  In  the 
table  B2D(I)  1=1,36. 

These  move/modifications  are  shown  In  Figure  3-8.  The  three 
labels  RECADR,  LOC,  B2D  Identify  the  names  of  the  three  tables  by 
which  the  data  are  later  accessed.  All  the  remaining  disk  I/O  and 
array  parameters  are  now  generated  from  the  other  preprocessing 
parameters . 
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Figure  3-7.  HSDB  Preprocess  Parameters 


All  the  level  0,1, 2, 3 routines  for  each  processing  pass  are 
obtained  from  the  object  module  associated  with  that  process. 

For  example.  Figure  3-9  shows  the  listing  of  the  object  module 
associated  with  Area  Fill  Processing.  PASSCT  defines  the  number 
of  complete  processing  passes  required  to  generate  the  Area  Fill. 
LEVnm  defines  the  tables  of  the  level  m routines  within  the 
particular  processing  pass  number  n^. 

The  label  names  associated  with  each  of  these  levels  define 
the  order  of  processing,  and  the  name  of  each  routine  within  each 
level.  For  example,  the  first  (and  only)  Level  0 routine  to  be 
executed  is  EDITSET. 

The  routines  at  Level  3 (first  pass)  are  the  Load  Array 
Routine  LDAl,  Line  Thinning  Routine  THIN2  etc. 

At  the  3rd  Level  the  appended  number  1 or  2 identifies  which 
Page  the  specific  routine  will  be  executed  from.  Notice  that  this 
is  an  alternating  process  between  Page  1 and  Page  2. 

The  two  specific  labels  NXTLVL,  and  X'80000000'  have  partic- 
ular meanings.  NXTLVL  tells  the  executive  to  go  to  the  next 
1 ower  level  routines  and  proceed  until  another  NXTLVL  is  found  or 
the  label  X'80000000'  is  located.  X'80000000'  tells  the  executive 
to  return  to  the  next  higher  level  routine  and  proceed  by  checking 
the  name  following  the  NXTLVL. 

Any  working  tables  and  buffer  areas  that  may  be  required 
during  processing  are  defined  as  subsequent  labels  within  the  module. 

The  level  0,1,  and  2 routines  are  processed  similarily  by 
the  Executive.  However,  the  processing  of  the  Level  3 routines 
requires  further  description. 

The  Level  3 routines  are  executed  in  order  by  alternating 
their  loading/execution  from  Page  1 and  Page  2,  Page  0 is  reserved 
for  the  execution  of  arithmetic  subroutines. 
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Figure  3-9.  Area  Fill  Processing  Module 
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The  actual  operations  involved  in  executing  the  Level  3 
routines  are  as  follows: 

1.  Initialize  SUBPASS3. 

2.  Set  the  Array  Input/Output  buffer  addresses  RESETLOC(I), 
1=1,28  in  a spare  table  LOC(I),  1=1,28. 

3.  Load  the  first  Level  3 routine  into  Page  1. 

4.  Initiate  loading  of  the  next  routine  into  the  second  Page 
(Page  2 initially),  and  execute  the  other  routine. 

5.  After  the  routine  is  executed,  repeat  step  4. 

6.  Repeat  steps  4 and  5 until  the  marker  X'80000000'  (end 
of  Level  3 routines)  is  found. 

7.  Execute  the  last  routine. 

8.  Decrement  SLIBPASS3. 

9.  Repeat  steps  2 through  8 until  SUBPASS3  is  zero. 

There  is  another  condition  that  prevents  the  Level  3 and 
Level  2 routines  from  exiting  to  the  next  higher  level.  The 
Level  2 routines  Load/Unload  a data  buffer  in  Bulk  Core  Memory. 

The  size  of  this  area  is  determined  by  the  amount  of  space  required 
for  tables  etc.,  during  processing.  This  area  is  equivalent  to  N 
array  loads  of  data  where: 

N array  loads  = 2 x number  of  4K  buffer  areas  allocated 

Therefore,  for  an  N array  load  of  data  residing  in  buffer  memory 
the  number  of  subpasses  at  level  3 (SUBPASS3)  to  completely  process 
this  data  will  be  jParray s^ i'n  "system » where  both  of  these  variables 
is  some  multiple  of  2.  Similarily  the  number  of  subpasses  at 
level  2 (SUBPASS2)  will  be: 

_ # vertical  scan  spots/record  . 

# vertical  scan  spots  in  N-arfay  load 

This  parameter  was  previously  generated  in  the  Preprocessing  routine 
and  is  equivalent  to  the  number  of  reformat  loads/buffer  load. 
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3.4  TAPE/DISK  ROUTINES 
3.4.1  Function  ’ 

The  function  of  the  Tape/Disk  routines  are  twofold. 

i)  To  control  the  Tape  Input/Output  routines. 

ii)  To  reformat  the  unprocessed/processed  data  and  trans- 
fer the  reformatted  data  between  the  reformat  area  in 


Bulk  Core  and  the  RK05  disk. 


3.4.2 


The  approach  was  to  write  two  routines  to  perform  the  two 
way  data  transfers  and  Input/Output  requests. 

3.4.3  Impi ementati on 

The  two  routines  are  written  in  APPLE  and  are  executed  as 
Level  1 routines  out  of  Bulk  Core.  The  reformatting  subroutines  are 
executed  out  of  page  memory  to  reduce  execution  time. 

3.4.3. 1 General  - The  routine  that  calls  the  tape  input  routine, 
and  performs  that  input  reformatting  is  called  TPRDF.  The  tape 
output  routine  is  called  TPWF. 

3. 4. 3. 2 TPRDF  - This  routine  controls  the  tape  input  routine  and 
executes  a subroutine  called  B2DI  which  reformats  the  input  records 
and  offloads  them  to  the  RK05  disk.  TPRDF  performs  the  following 
initial  operations: 

•Loads  the  B2DI  routine  into  Page  1 for  subsequent  execution. 
•Set  the  index  pointer  used  by  D2BI  DISCINDXI  to  zero. 


•Sets  the  initial  disk  block  number  at  4300. 


The  disk 


blocks  following  will  contain  the  overlap  records  which  will 
be  the  first  records  of  the  next  data  block  (BSDBLK). 

Sets  the  starting  disk  block  numbers  for  the  multiple  array 
loads/unloads  during  the  DISK-BUFFER  I/O  (BLKB2D). 
Initializes  the  number  of  plot  records  in  each  buffer  load 
(RECBF). 

Sets  the  number  of  overlap  records  between  each  data  block 
(OVLP). 


•Initializes  the  pointer  (SRECAD)  to  the  table  - (RECADR) 
which  contains  the  Bulk  Core  Memory  base  addresses  of  where 
the  Tape  Input  records  will  be  loaded. 

•Sets  the  number  of  unique  input  plot  records/data  block  load 
(IRCM4). 

•Sets  the  total  number  of  input  plot  records/data  block, 
including  the  overlap  records  (IPRBLK), 

•Sets  the  total  number  of  plot  format  records  expected  to  be 
input  or  generated  for  the  overlay  sheet  (NREC). 

•Sets  a variable  (EOFCK)  during  the  first  and  the  last  data 
blocks. 

EOFCK  = % of  number  of  overlap  records  for  first  block, 
EOFCK  = 300  for  last  block. 

For  the  first  data  block  EOFCK  is  used  as  a decrementabl e count  to 
generate  blank  records  in  the  first  half  of  the  overlap  area.  For 
the  last  data  block  EOFCK  is  used  as  a decrementabl e count  to 
generate  blank  records  to  fill  up  the  last  data  block  for  process- 
ing. 

The  basic  sequence  of  operations  in  TPRDF  are  as  follows; 

1.  Call  the  routine  to  read/generate  one  plot  tape  format 
record. 

2.  Repeat  step  1 until  buffer  is  full;  i,e.,  RECBUF=0. 

3.  Reformat  a section  (multiple  array  load)  of  the  tape 
records  via  reformat  area  to  the  disk. 

4.  Repeat  step  3 until  all  tape  records  completely  off- 
loaded to  the  disk. 

5.  Repeat  steps  1-4  until  all  input  records  for  a data 
block  have  been  read;  i.e.,  IRCM4=0. 

Steps  3 and  4 are  the  subroutine  B2DI  which  will  be  discussed 
subsequently . 
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Routine  TPRDF  also  executes  to  subroutines  which  save  the 
overlap  records  (SAVREC)  and  retrieves  them  at  the  start  of  the 
next  data  block  (GETREC).  A set  of  disk  blocks  are  reserved  for 
these  overlap  records,  the  initial  disk  block  being  defined  as 
4300jg. 

The  subroutine  B2DI  is  called  by  TPROF  and  reformats  the  tape 
input  data  onto  the  RK05  disk  for  subsequent  load  by  the  Disk-to- 
Buffer  Input  (D2BI)  routine  described  in  subsection  3,6,  The 
B2DI  subroutine  transfers  the  tape  buffer  data  to  disk  in  N-array 
loads  (N  is  a multiple  of  4).  This  subroutine  is  a Level  1 sub- 
routine of  TPRDF,  but  is  executed  out  of  Pages  to  increase  efficiency, 
The  steps  involved  in  B2DI  are  as  follows; 

1.  Perform  the  following  initialization. 

Generate  the  table  RECADRL  from  the  table  of  the  tape 
record  base  addresses  RECADR. 

Set  up  the  TRANBLOK  to  point  to  the  REFORMAT  AREA  and 
for  output  to  the  disk. 

Set  the  number  of  N-array  loads  (N8ARYL). 

Set  the  disk  block  index  pointer  to  zero  (DISCINDX). 
Initialize  the  record  address  table  index  to  zero. 

2.  Get  the  number  of  tape  records/buffer  load  (NRECBL). 

3.  Get  the  next  record  address  from  RECADRL. 

4.  Transfer  the  number  of  32-bit  words  corresponding  to 
N-array  loads  (of  one  record)  to  the  reformat  area, 

5.  Decrement  the  pointer  to  the  last  word  transferred  from 
the  tape  record  by  the  number  of  words  in  the  overlap 
area . 

6.  Save  this  pointer  in  RECADRL. 

7.  Repeat  steps  3 through  6 until  NRECBL=0. 

8.  Get  the  base  disk  block  number  from  BLKB2D  (DP)  where  DP 
is  the  index  pointer  DISCINDX.  Set  the  disk  block  number 
in  TRANBLOK. 
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9.  Output  the  Reformat  Area  to  the  disk, 

10.  Update  the  base  disk  block  number  by  the  number  of  disk 
blocks  just  transferred  (NBLKTRN). 

11.  Update  the  base  disk  block  table  index. 

12.  Repeat  steps  2 through  11  until  all  N-array  loads  have 
been  reformatted  and  output  to  the  disk  (N8ARYL=0). 

13.  Return  to  TPRDF. 

3. 4. 3. 3 TPWF  - This  routine  controls  the  tape  output  routine  and 
executes  a subroutine  called  D2B0  which  reformats  the  processed  data 
into  output  (plot  format)  tape  records  from  the  RK05  disk. 

TPWF  performs  the  following  initialization  operations; 

• Sets  the  index  pointer  DISCINDXO  to  zero.  This  pointer  is 
used  by  B2D0  to  index  the  buffer  to  disk  table  which  contains 
the  base  disk  block  numbers  which  will  contain  the  output 
data  for  reformatting  (B2D), 

• Initiates  the  movement  of  the  Di sk-to-Buf fer  Output  routine  to 
Pages . 

• Sets  the  table  B2D  in  secondary  storage  defined  as  BLKB2D. 

• Sets  the  number  of  output  plot  records/data  block  (RECBLK). 

TPWF  then  performs  the  following  steps: 

1.  Sets  the  number  of  output  records/buffer  load, 

2.  Sets  the  index  pointer  SRECAD  to  the  table  of  output 
record  base  addresses  in  RECADR. 

3.  Performs  the  subroutine  D2B0. 

4.  Calls  the  routine  to  generate/write  one  the  equivalent 
of  one  plot  format  record. 

5.  Performs  step  4 until  records/buffer  is  zero. 

6.  Repeats  steps  1 through  5 until  RECBLK  is  zero. 

7.  Exits  Routine. 

The  subroutine  D2B0  is  called  by  TPWF  and  reformats  the 
output  processed  data  from  the  Buf fer-to-Di sk-Output  routine  B2D0 


44 


rntmi  ra  I 


described  in  section  3.6.  The  D2B0  subroutine  transfers  the 
I processed  data  from  the  disk  to  the  reformat  area  in  N-array 

I loads.  This  subroutine  as  a Level  1 subroutine  of  TPWF  which  is 

executed  out  of  Pages  to  increase  efficiency.  The  steps  involved 
I in  D2B0  are  as  follows: 

. 1.  Perform  the  following  initialization. 

I • Generate  the  table  RECARDRl  from  the  table  of  the  tape 

record  base  addresses  RECADR. 

I • Set  up  the  TRANBLOK  to  point  to  the  REFORMAT  Area  and 

for  input  from  the  disk. 

T • Set  the  disk  block  index  pointer  to  zero  (DISCINDX2). 

• Set  the  number  of  N-array  loads  (N8ARYL). 

I 2.  Set  the  number  of  plot  tape  records/buffer  load  (NRECBL). 

• 3.  Set  the  initial  disk  block  number  in  the  TRANBLOK  for 

I the  disk  output  from  the  disk  to  buffer  block  table  BLKD2B, 

1 using  the  index  pointer  DISCINDX2. 

4.  Transfer  from  disk  to  reformat  area  the  disk  blocks 

I corresponding  to  the  processed  data  to  be  reformatted 

for  output. 

5.  Update  the  BLKD2B  disk  block  number. 

""  6.  Reformat  the  data  into  the  partial  tape  record. 

7.  Update  the  start  address  of  the  RECADRl  table  for  the 

••  next  reformatting. 

8.  Repeat  steps  5 through  7 until  NRECBL  is  zero. 

9.  Repeat  steps  2 through  8 until  all  of  the  N-array  loads 
have  been  output  and  the  tape  records  are  generated. 

10.  Return  to  TPWF. 

3.5  TAPE  I/O  ROUTINES 
i • 3.5.1  Function 

The  function  of  these  routines  is  to  perform  the  transfer  of 
the  input  data  and  the  processed  output  data  between  the  tape 
devices  and  the  Bulk  Core  buffer  area. 


3.5.2  Approach 

The  approach  was  to  write  several  APPLE  routines  to  allow 
the  input/output  of  raster  data  in  several  formats.  These  routines 
are  linked  directly  to  the  Tape/Disk  routines.  The  different  input 
formats  are: 

1.  E.T.L  run- 1 ength-coded  scanned  data. 

2.  E.T.L  binary  plot  format  data. 

3.  D.M.A  run- 1 ength-coded  scanned  data. 

The  output  formats  are: 

1.  E.T.L  binary  plot  format  data. 

2.  D.M.A  run-1 enght-coded  plot  data. 

3.5.3  Impl ementati on 

The  Tape  I/O  routines  are  written  in  APPLE  and  executed 
at  Level  1 as  subroutines  of  TPRDF,  and  TPWF,  the  Tape/Disk  routines. 

3. 5. 3.1  General  - The  Tape  I/O  routines  available  are  the 
fol 1 owi ng  : 

• ETLIOl  - This  routine  reads  ETL  binary-plot  format-data, 

and  writes  ETL  binary-plot  format  data. 

• ETLI02  - This  routine  reads  ETL  run- 1 ength-coded  format 

data,  and  writes  ETL  binary-plot  format  data. 

• DMAIO  - This  routine  reads  DMA  run- 1 ength-coded  format 

data,  and  writes  DMA  run-length-coded  format  data. 

• DMAETL  - This  routine  reads  DMA  run- 1 ength-coded  format 

data,  and  writes  ETL  binary-plot  format  data. 

Any  one  of  these  I/O  subroutines  may  be  selected  prior  to  execution 
of  the  preprocessing  routine  by  the  selection  on  the  lOTYPE  para- 
meters card  discussed  in  subsection  3.2.  After  execution  of  the 
preprocessing  routines  will  come  a command  to  load  the  selection 
I/O  routine  prior  to  load  and  executing  the  particuliar  processing 
routines.  Each  of  the  I/O  routines  is  loaded  at  location  8200^g  in 
Bulk  Core  Memory.  The  first  instruction  check  as  switch  (Data 
Pointer  Register)  which  is  set  by  the  calling  Tape/Disk  routine, 
to  select  either  the  input  or  the  output  part  of  the  process. 

Figure  3- 10  shows  the  top-level  flow  of  the  I/O  routines. 
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Each  of  the  different  Tape  I/O  routines  is  discussed  in  detail 
in  the  following  pages. 

3. 5. 3. 2 ETLICl  - This  routine  reads  ETL  binary-plot  formatted 
data,  and  writes  ETL  binary-plot  formatted  data.  Fi gure  3- 1 1 shows 
the  flow  diagram  for  this  routine.  This  routine  uses  the  direct 
transfer  functions  used  in  the  Tape/Disk  routines.  The  TRANBLOK 
instructions  are  set  up  prior  to  execution  of  the  data  transfer. 

The  buffer  area  that  will  contain  the  input  record  is  cleared 
ahead  of  the  actual  transfer.  Completion  of  the  input  or  output 
transfer  results  in  a return  to  the  calling  routine. 

3. 5. 3. 3 ETLI02  - This  routine  reads  ETL  run- 1 ength-coded  (RLC) 
data,  and  writes  ETL  binary-plot  formatted  data.  The  output 
portion  of  this  routine  is  identical  to  that  discussed  in  the 
previous  subsection.  The  input  portion  of  this  routine  reads  in 
variable  length  run- 1 ength-coded  records  and  reformats  each  group 
of  eight  records  into  the  binary  format  equivalent.  Control  is 
then  returned  to  the  calling  routine. 

For  an  explanation  of  the  ETL  run- 1 ength-coded  format  the 
reader  is  referenced  to  Appendix  A of  this  report.  Figure  3-12 
shows  the  top-level  flow  diagram  for  inputting  ETL  RLC  data. 

3. 5. 3. 3.1  General  - Basically  eight  RLC  input  records  (scan  lines) 
are  required  to  generate  one  binary  plot  record.  These  records 
are  alternately  input  to  one  of  two  buffers.  The  first  buffer 
starts  at  location  FOOOjg,  and  the  second  buffer  starts  at  loca- 
tion FBOOjg.  Each  buffer  area  is  2048  words.  After  the  first  of 
the  eight  records  has  been  input  the  reformatting  of  each  record 
into  the  binary  plot  format  is  performed  while  the  next  record  is 
being  transferred.  This  'double  buffering'  procedure  absorbs 
much  of  the  reformatting  time. 

After  each  set  of  eight  records  has  been  input  and  reformatted 
to  generate  one  plot  record,  control  is  returned  to  the  calling 
routine  (TPRDF). 
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Figure  3-11.  ETLIOl  Flow  Diagram 
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Figure  3-12.  ETL  R.L.C.  Input  Top-Level  Flow  Diagram 


3.5. 3. 3. 2  Reformatting  - Figure  3-13shows  the  flow  diagram  for 
the  reformatting  of  the  RLC  data  to  binary  plot  data.  Essentially, 
the  transitions  from  OFF  to  ON  (Wordl,  Bytel),  and  from  ON  to  OFF 
(Word2,  Byte2)  are  computed.  The  bit  within  Byte2  of  the  second 
transition  word  M0RD2  is  then  set.  This  bit  depends  upon  the  scan 
number  (0,1,2. ..7)  currently  being  input.  The  same  bit  within  all 
bytes  starting  with  WORDl,  BYTEl  are  then  set  in  a sequential 
manner  until  this  same  bit  is  located.  The  process  is  then  repeated 
for  all  pairs  of  transitions. 

3. 5. 3. 4 DMAIO 

3. 5. 3. 4.1  Function  - The  function  of  the  DMA  I/O  routines  is  to 
convert  DMA  run- 1 ength-coded  data  to  binary  for  processing  and  to 
convert  processed  binary  data  back  to  run-length. 

3. 5. 3. 4. 2 Approach  - Because  I/O  is  an  area  most  likely  to  change 
in  a production  system,  the  approach  was  to  convert  the  data  to 
and  from  binary  plot  tape  format  and  then  utilize  existing  tape  and 
disk  I/O  routines. 

3. 5. 3. 4. 3 Implementation  - There  are  two  basic  routines  involved 
in  the  conversion  of  DMA  run  length  coded  (RLC)  data.  DMARLI 
converts  RLC  data  to  binary  data.  DMARLO  converts  binary  data  to 
RLC  data.  The  format  of  the  DMA-RLC  data  is  given  in  Appendix  A. 

The  implementation  of  these  routines  are  discussed  next. 

The  salient  tasks  of  the  SMARLI  routine  are  outlined  below: 

1.  Transfer  (TRAN)  in  1 RLC  record.  (Eight  RLC  records 
are  required  to  create  one  binary  record.) 

2.  Using  the  HI  WC,  LO  WC  pair,  load  the  field,  RLC, 
with  the  run-length  values.  Load  until  all  values 
have  been  loaded,  or  until  there  is  no  room  left  in 
the  array. 

3.  Call  VLSUM$  to  calculate  running  sums  for  each  array. 

4.  Add  last  sum  of  each  array  to  all  sum  values  in 
succeeding  arrays. 
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Figure  3-13.  Reformatting  (R.L.C.  to  Binary  Plot) 
plow  Diagram  (Sheet  2 of  2) 


5. 


Lastly,  add  SAVSUM  (any  sum  value  from  prior  array 
loads  this  record)  to  the  sum  values. 

6.  The  data  in  the  SUM  field  now  contains  the  cumulative 
distance  from  the  start  of  the  sheet  to  the  transition 
flagged  by  the  most  significant  bit  of  RLC. 

7.  Calculate  the  following  parameters: 

START  (offset  in  Binary  record)  = SUM-RLC 
start  BYTE  = Remainder  (START/4) 

Word  address  = value  pointed  to  by  SRECAD 
+ Quotient  (START/4) 

NSLTS  (number  of  slots  available)  = 4-BYTE 
For  all  RLC  > NSLTS: 

FWC  (first  word  count)  = NSLTS 
WC  (full  word  count)  = Quotient  ( RLC-NSLTS )/4 
LWC  (last  word  count)  = Remainder  ( RLC-NSLTS )/4 
For  all  RLC  < NSLTS 
FWC  = RLC 
WC  = 0 
LWC  = 0 

8.  The  bit  to  be  set  in  the  word  is  determined  by  the 
least  significant  three  bits  of  the  run-length  record 
address  (HI  Y,  LO  Y).  These  three  bits  index  an  8 
word  table,  BITTAB  (refer  to  Figure  3-14).  The  proper 
pattern  is  selected  and  stored  in  BITPAT. 

If  record  number  = ^^10 
(Binary  = 011010) 

Then  least  significant  3 bits 
= 010 

BITPAT  = X'20202020' 


X'80808080' 
X' 40404040' 
X'20202020' 
X' 10101010' 
X' 08080808' 
X'04040404' 
X'02020202' 
X'OIOIOIOI' 


Figure  3-14.  BITTAB  format 
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9.  Once  steps  7 and  8 have  been  completed,  the  proper 
bits  are  set  in  the  binary  record  by  calling  SETNBITS 
and  SETMORDS  (described  below). 

10.  When  all  bits  have  been  set  this  arrayload,  DONE 
flag  is  checked.  If  DONE^l,  steps  2-9  are  repeated 
until  D0NE=1. 

11.  When  D0NE=1,  the  COUNT  is  checked  to  see  if  all  8 
RLC  records  have  been  incorporated  into  the  binary 
record.  Steps  1-11  are  repeated  until  C0UNT=0. 

12.  The  program  then  exits. 

The  SETNBITS  routine  is  called  when  FWC  or  LWC  are  greater 
than  zero.  It  sets  the  number  of  bits  in  a word  as  specified  in 
the  BL  (FWC  or  LWC),  starting  at  the  byte  pointed  to  by  the  DP 
(BYTE).  The  bytes  of  the  F-register  are  initially  all  set.  This 
is  achieved  by  loading  F with  BITPAT.  Depending  upon  the  start 
byte  and  the  number  of  bytes  to  be  set,  the  corresponding  8-bit 
registers  of  the  F-register  are  cleared.  The  contents  of  the 
F-register  is  then  .OR.ed  with  the  word  from  the  binary  record 
table  (pointed  to  by  ADR)  and  the  resultant  word  stored  back  into 
the  binary  record  table.  ADR  is  incremented  prior  to  the  return. 

The  SETWORDS  routine  is  called  whenever  WC  is  greater  than 
zero.  It  gets  WC  words  from  the  binary  table,  .OR.s  them  with 
BITPAT,  and  stores  the  results  back  into  the  binary  table.  ADR  is 
bumped  by  WC  prior  to  the  return. 

The  salient  tasks  of  the  DMARLO  routine  are  outlined  below: 

1.  Upon  entry,  the  program  performs  certain  initialization 
tasks  that  must  be  done  for  each  binary  record.  It 
also  creates  some  flag  bits  which  are  required  later 
on  in  the  conversion  process.  The  start  address  of 
the  binary  record  and  its  length  are  stored  in  the 
location  LNGADR. 
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2.  Data  are  then  loaded  from  bulk  core  bytewise  into 
the  array  field  BIN. 

3.  Eight  run-length  records  are  created  from  one 
binary.  To  avoid  reloading  the  binary  data  8 times, 
the  run-length  records  are  built  and  maintained  in 
the  array  in  fields  RECn,  n=l-8.  To  simplify  the 
programming  however,  all  run-length  records  are 
created  and  output  from  the  same  field,  BLDAOPT. 

Prior  to  processing,  the  appropriate  RECn  field  is 
moved  to  BLDAOPT.  This  field  is  moved  back  after 
process i ng . 

4.  Transition  points  are  found  by  shifting  the  binary 
slice  (Bits  0-7  of  BIN)  up  one  and  exclusive  .OR.ing 
the  data  with  itself. 

5.  The  black/white  flag,  RLC,  is  created  by  AND.ing  the 
transition  flag  with  the  binary  source  slice.  The 
run-length  value  (RIVAL)  is  then  equal  to  the  address 
of  the  current  transition  minus  the  address  of  the 
previous  transition  value. 

6.  The  BLDREC  is  called  to  build  a run-length  record. 

7.  Steps  5 and  6 are  repeated  for  all  transition  points 
in  a slice. 

8.  Steps  3 thru  7 are  repeated  for  each  run-length 
record  in  the  binary  record. 

9.  Steps  2 thru  8 are  repeated  until  the  entire  binary 
record  has  been  converted.  The  subroutine  CLEANUP  is 
then  called  and  the  program  returns. 

The  BLDREC  adds  the  data  being  created  in  RLVAL  to  the  record 
data.  The  set-up  for  the  call  to  BLDREC  is  shown  below. 

CH=Quotient  (RLVAL/128)  with  MSB=RLC 
CL=Remainder  (RLVAL/128)  with  MSB=RLC 
BL=CH  without  RLC 
DP=CL  without  RLC 
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If  BL=0,  RIVAL  is  less  than  or  equal  to  127  so  only  the  contents 
of  CL  need  be  added  to  BLDAOPT.  If  BL^O,  the  contents  of  CH  must 
be  added  to  the  record  along  with  BL-number  of  7F  or  FF  values. 

If  DP=0  (i.e.,  RLVAL  is  an  exact  multiple  of  128),  CL  need  not  be 
stored . 

BLDREC  must  check  prior  to  each  store  to  insure  that  there 
is  room  in  BLDAOPT.  If  there  is  no  room,  OPTREC  is  called  to 
output  BLDAOPT  to  the  proper  record  buffer. 

The  OPTREC  unloads  the  BLDAOPT  field  into  a buffer  in  bulk 
core.  Using  BIN  as  a pointer,  the  address  of  the  record  is 
retrieved  from  RECPTR  and  the  half  word  count  is  retrieved  from 
BYTCT. 

The  OPTREC  outputs  BLDAOPT  and  counts  the  number  of  half  words 
being  stored.  The  common  register  is  originally  padded  with  7F 
values  to  insure  that  an  even  number  of  run  length  values  are 
stored  away. 

Once  the  data  have  been  stored,  RECPTR  and  BYTCT  are  returned 
to  memory.  BLDAOPT  is  cleared  and  the  program  returns. 

The  CLEANUP  first  insures  that  all  data  have  been  output, 

I if  not,  OPTREC  is  called.  The  header  data  (HI  WC,  LO  WC)  (HI  Y, 

: LO  Y)  are  then  added  to  each  record  and  the  records  are  .TRAN.ed 

to  tape. 

■ 

3.6  DISK  I/O  ROUTINES 
3.6.1  Functi on 

The  disk  input/output  routines  transfer  data  between  the 
RK05  disk  and  a Bulk  Core  Memory  buffer. 

r 

; 3.6.2  Approach 

I The  approach  was  to  write  two  routines  in  APPLE  to  perform 

the  two  way  transfers. 

3.6.3  Impl ementati on 

I The  two  routines  are  written  in  APPLE  and  are  executed  as 

Level  2 routines  from  Buld  Core  Memory. 
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3. 6. 3.1  General  - The  two  routines  transfer  a block  of  arrays  of 
data  between  the  buffer  and  the  disk.  The  number  of  arrays  of 
data  depends  on  the  size  of  the  buffer  area  being  used,  and  is 
normally  a multiple  of  four.  The  routine  that  transfers  the  un- 
processed data  from  the  disk  to  the  buffer  is  called  D2BI.  The 
routine  that  transfers  the  processed  data  from  the  buffer  to  the 
disk  is  called  B200. 

3.6. 3.2  Pi sk-to-Buffer  Transfer  - This  routine  transfers  reformatted 
unprocessed  data  generated  by  the  Tape/Disk  Routines  to  a buffer 

in  Bulk  Core.  The  D2BI  routine  performs  the  following  operation: 

1.  Sets  up  the  TRANBLOK  parameters  associated  with  the  data 
transfer,  namely, 

i)  Buffer  Location, 

ii)  Number  of  words  to  transfer, 

iii)  Direction  of  transfer, 

iv)  The  starting  disk  sector  number  of  the  data  block 
that  is  to  be  transferred 

2.  Performs  the  di s k-to-buf f er  transfer,  and 

3.  Exits  the  routine  after  completing  the  data  transfer. 

3. 6. 3. 3 Buf f er- to-Di s k Transfer  - This  routine  transfers  processed 
data  output  from  the  arrays  and  located  in  the  data  buffer  used 

by  D2BI  to  the  disk.  The  B2D0  routine  performs  the  following 
operati ons : 

1.  Sets  up  the  TRANBLOK  parameters  associated  with  the  data 
transfer,  namely, 

i)  Buffer  location, 

ii)  Number  of  words  to  transfer, 

iii)  Direction  of  transfer, 

iv)  The  starting  disk  sector  number  for  the  data 
block  being  output, 

2.  Performs  the  buff er-to-di sk  transfer,  and 

3.  Exits  the  routine  after  completing  the  data  transfer. 
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3.7  ARRAY  LOAD/UNLOAD  ROUTINES 

3.7.1  Function 


1 


The  array  load/unload  routines  transfers  the  preprocessed/ 
processed  raster  data  between  Bulk  Core  Memory  and  the  STARAN  arrays. 

3.7.2  Approach 

The  approach  was  to  reunite  those  routines  already  written 
for  Loading/unloading  the  arrays  to  take  advantage  of  mix-mode 
operations  and  thus  speed  up  this  process. 

3.7.3  Implementati on 

The  modified  array  input/output  routines  are  written  in 
APPLE  assembly  language  and  machine  code,  and  are  executed  as 
Level  3 routines  from  the  Page  memories. 

3. 7. 3.1  General  - These  two  routines  load/unload  the  arrays  from/to 
Bulk  Core  Memory.  These  routines  utilize  three  subroutines, 

LDBYT$,  ORDERS,  and  STBYT$  which  were  developed  to  efficiently 
transfer  byte  values  between  control  memory  and  array  modules. 

These  subroutines  are  described  in  subsection  3. 7. 3,4. 

3. 7. 3. 2 Load  Arrays  - The  load  array  routine  loads  the  arrays 

with  raster  data.  This  data  has  been  extracted  from  the  reformatted 
data  stored  on  the  RK05  disk  using  the  Di s k-to-Buffer- Input  (D2BI) 
Routine  described  in  subsection  3.6.  The  operations  required  to 
load  the  arrays  are  as  follows: 

1.  Perform  initialization  of; 

i)  Number  of  arrays  in  system, 

ii)  Number  of  8-bit  (byte)  array  fields  to  be  loaded, 

2.  Get  pointer  to  Bulk  Core  Memory  containing  the  64  words 
to  be  stored  into  the  arrays, 

3.  Get  array  field  pointer. 

4.  Load  the  array  field  using  LDBTY$* 
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4.  Load  the  array  field  using  LDBYT$*, 

5.  Repeat  step  4 for  all  arrays, 

6.  Save  Bulk  Core  Memory  pointer  for  next  time, 

7.  Repeat  steps  2 through  6 for  all  array  fields, 

8.  Order  the  flipped  data  using  ORDERS,  and 

9.  Exit  Routine. 

*It  should  be  noted  that  each  32-bit  (4  byte)  word  in  Bulk 
Core  is  stored  in  reverse  order.  Therefore,  each  four  bytes 
of  data  is  flipped  to  its  correct  location  within  each  32-bit 
section  of  the  X-response  store  prior  to  loading  into  the 
arrays  as  shown  in  Figure  3-15. 

3. 7. 3. 3 Unload  Arrays  - The  unload  array  routine  unloads  the 
processed  array  data.  This  data  will  later  be  stored  on  the  RK05 
disk  using  the  Buffer-to-Disk- Output  Routine  {B2D0)  described  in 
subsection  3.6.  The  operations  required  to  unload  the  arrays  are 
as  follows: 

1.  Perform  initialization  of; 

i)  Number  of  arrays  in  system, 

ii)  start  column  of  data  to  0/P, 

iii)  Number  of  words/8-bit  field  to  output, 

2.  Move  output  data  area  to  the  top  left  edges  of  the  arrays, 

3.  Determine  number  of  output  words/8-bit  field, 

4.  Get  next  field  pointer  and  Bulk  Core  Memory  output  area 
pointer , 

5.  Unload  one  array  field  of  data  from  all  arrays  to 
Bulk  Core  Memory, 

6.  Save  Bulk  Core  Memory  output  area  pointer, 

7.  Repeat  steps  4 thru  6 until  all  8-bit  fields  have  been 
output,  and 

8.  Exit  Routine. 

3. 7. 3. 4 LDBYT$.  ORDERS,  STBYT$  - These  three  subroutines  perform 
efficient  transfers  of  byte  values  between  Bulk  Core  memory  ad 
and  the  array  modules.  The  area  in  memory  is  a block  of  64  words. 
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Figure  3-15.  Data  Correction  During  Array  Loading 
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each  containing  four  bytes  or  a total  of  256  byte  values.  The 
most  significant  bit  address  of  the  array  field  must  be  a multi- 
ple of  8. 

3. 7. 3. 4.1  Load  Array  - The  load  array  (LDBYT$)  subroutine  will 
load  an  8 bit  field  of  a selected  array  module  with  the  256  byte 
values  contained  in  a 64  word  block  of  control  memory.  The  bytes 
will  be  stored  in  mixed  mode  order.  The  calling  sequence  is: 

FPl  = array  module  number 

FP2  = MSB  of  destination  field  (must  be  a multiple  of  8) 

DP  = address  of  1st  of  the  64  source  control  memory  words 
RO  = return  address 

The  following  registers  are  used  by  LDBYT$: 
X,C,FL1,FP1,FP2,BL,DP,ASL,R0 

The  execution  time  for  LDBYT$  is  120  microseconds  and  the  listing 
is  shown  in  Figure  3-16. 

Right  after  calling  LDBYT$,  the  routine  0RDER$,  discussed  next, 
should  be  called  to  put  the  byte  data  in  array  memory  into  the 
natural  order. 

3. 7. 3. 4. 2 Order  Data  - The  Order  Data  (ORDERS)  subroutine  will 
transform  the  array  memory  byte  data  from  mixed  mode  order  into 
natural  order  for  subsequent  array  operations.  If  the  array 
memory  byte  data  is  already  in  natural  order,  ORDERS  will  transform 
the  data  into  mixed  mode  order  so  STBYT$  can  properly  transfer 

the  data  from  array  memory  into  control  memory.  The  calling 
sequence  is: 

AS  = array(s)  to  be  modified 
FL2  = length-1  of  the  field 
FPE  = MSB  position  of  the  field 
RO  = return  address 

The  following  registers  are  used  by  ORDERS: 
X,M,C,FL1,FPE,FL1,FPE,FP1,R0,R1 
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ASL.OP 
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0002 

0002 

38014807 

GEN. 32 

X*38014807* 

PUT  AMRl  IN  BYTE  MODE 

0003 
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3F030014 
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LOOP. 4 
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36020040 
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0005 
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0006 
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0007 

0007 
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X(l) 

0008 

0008 
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LR 

C.16(DP) 

0009 

0009 
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X(2> 

OOOA 

OOOA 

36020018 

LR 

Ct24(0P) 

OOOB 

OOOB 
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X(3I 

OOOC 

OOOC 

36020020 
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C.32(0P) 

OOOD 

OOOD 

01C499A0 

sew 

X(4) 

OOOE 

OOOE 

36020028 

LR 

C.40(OP) 

OOOF 

OOOF 

01C599A0 

sew 

X(5) 

0010 

0010 

36040030 

LR 

C«48(0P)«2 

INCR  OP 

0011 

0011 

01C699A0 

sew 

X(6) 

0012 

0012 

36050037 

LR 

C«55(0P).3 

iNCR  OP.DECR  BL 

0013 

0013 

41C799A0 

sew 

X(7) 

0014 

0014 

1FE00003 

LDBYTS2 

s.w 

X.FP12^ 

STORE  XlINCE  FP2 

0015 

0015 

3040000A 

LRR 

DP.ASL 

RESTORE  OP 

0016 

0016 

29010019 

BZ.BL 

LDBYTS3 

EXIT  IF  DONE 

0017 

0017 

28140001 

INCR 

OP 

ELSE  INCR  DP 

0018 

0018 

28010003 

B 

LDBYTSI 

AND  REPEAT  OUTER  LOOP 

0019 

0019 

38014eFF 

LDBYTS3 

GEN. 32 

X'380148FF* 

RESTORE  AMRl  TO  WORC 

OOlA 

OOlA 

28080000 

B 

0(R0) 

RETURN  TO  CALLER 

f 


Figure  3-16.  Load  Array  (LDBYT$)  Listing 
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The  execution  time  for  0RDER$  is  1.25*  (field  length-1)  + 13 
microseconds  and  the  listing  is  shown  in  Figure  3-17. 

3. 7. 3. 4. 3 Dump  Array  - The  Dump  Array  (STBYT$)  subroutine  will 
transfer  an  8 bit  field  of  array  data  to  a block  of  64  words  of 
control  memory.  In  other  words,  the  opposite  function  is  perfor- 
med by  LDBYT$.  The  address  of  the  MSB  of  the  data  to  be  transferred 
must  be  a multiple  of  eight  and  the  data  must  be  in  mixed  mode 
order.  The  calling  sequence  is: 

FPl  = array  module  number 

FP2  = MSB  of  source  field  (must  be  a multiple  of  8) 

DP  = address  of  the  1st  of  the  64  destination  control 
memory  words 
RO  = return  address 

The  following  registers  are  used  by  STBYT$: 
X,C,FL1,FP1.FP2,BL,DP,ASL,R0 

The  execution  time  for  STBYT$  is  134  microseconds  and  the  listing 
is  shown  in  Figure  3-18. 

3.8  LINE  THINNING  ROUTINE 

3.8.1  Functi on 

The  line  thinning  routine  is  an  AAP  subroutine  which  thins 
line  data  to  a unit  (one-all)  thick  line. 

3.8.2  Approach 

The  approach  was  to  utilize  the  existing  routine,  modified 
to  automatically  terminate  thinning. 

3.8.3  Implementati on 

The  line  thinning  routine  (THIN)  is  a Level  3 routine 
written  in  APPLE  and  executed  in  Page  Memory. 

3.8. 3.1  General  - The  line  thinning  routine  was  modified  to 
automatically  terminate  when  all  lines  in  the  arrays  were  completely 
thinned  (i.e.,  no  more  calls  could  be  removed). 
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orders  start 

I 

(CHANGES  UORO  ORDER  OF  A FIELD  IN  ALL  ENABLED  ARRAYS 
(BETWEEN  NATURAL  ORDER  AND  ORDER  REQUIRED  BY  'LOBYTS' 
(AND  ‘STBYTS*  ROUTINES, 

I 

(USE  'ORDERS*  AFTER  'LDByTS'  TO  PUT  INPUT  BYTES  IN 
(NATURAL  ORDER. 


t 

(USE  'ORDERS'  BEFORE  'STBYTS'  SO  OUTPUT  BYTES  WILL 
(BE  IN  NATURAL  ORDER. 


(ENTRY;  AS=  ARRAYS  TO  BE  MODIFIED 

( FL2=  L'(FIEL0» 

( FPE=  M'JFIELO) 

I BAL«R0  ORDERS 

I 

(DISTURBS  X.M.C«FLltFL3»FPl.Rl 


t 

(TlMlN6:i.25(L'<FIELD)  )-fl3  MICROSEC 


t 

ENTRY 

ORDERS 

0000 

0000 

742055AA 

ORDERS 

LI 

CH.X'55AA' 

C=MASK  1 

0001 

0001 

720055AA' 

LI 

CL.X'55AA' 

0002 

0002 

33900009 

LI 

FP1.X»09' 

FLIP  CODE  (00001001) 

0003 

0003 

2CH0009 

BAL.Rl 

ORDERl 

call  inner  routine 

0004 

0004 

74203333 

LI 

CH.X'3333' 

C=MASK2 

0005 

0005 

7200CCCC 

LI 

CL.X'CCCC 

0006 

0006 

33900012 

LI 

FP1.X'12' 

FLIP  CODE  (OOOIOOICI 

0007 

0007 

2C110009 

BAL.Rl 

ORDERl 

CALL  INNER  ROUTINE 

0006 

0006 

26060000 

B 

0(R0) 

RETURN 

1 

(INNER 

ROUTINE  - SETS  UP  MASK  AND 

FLIPS  FIELD  between 

(MASKED 

WORDS. 

0009 

0009 

420086A0 

ORDERl 

SC.X'6eA0'  X(0) 

MASK  FROM  C 

OOOA 

OOOA 

40E099e3 

GENS .32 

X*40E099B3» 

REPLICATE  INTO  ALL  X 

OOOB 

OOOB 

40C099R3 

GENS. 32 

X'40C099B3' 

OOOC 

OOOC 

408099P3 

GENS. 32 

X'40a099B3' 

0000 

OOOD 

08000003 

L 

M.X 

M=MASK 

OOOE 

OOOE 

31000017 

LRR 

«FL1.FP3) . JFl2 

.FPE)  FIELD  PARAMETERS 

OOOF 

OOOF 

30000013 

LOOP 

0RDER2 

0010 

0010 

438008A5 

L 

X.FP3 

LOAD  X WITH  A BIT 

0011 

0011 

430000A3 

gens. 32 

X'430088A3' 

FLIP  X USING  FPl 

0012 

0012 

4B800001 

GENS. 32 

X'4B800001' 

SETUP  M 

0013 

0013 

53A00003 

0RDER2 

GENS. 32 

X'53A00003' 

STORE  X MASKE0{FP34 

0014 

0014 

00000001 

NOP 

0015 

0015 

28090000 

B 

0(R1) 

RETURN  TO  OUTER  ROUTE 

Figure  3-17.  Order  Data  (ORDERS)  Listing 


stbyts  start 

I 

ISTORCS  256  BYTES  FROM  AN  6>BIT  FIELD  OF  ALL  WORDS  IN  ONE 
(ARRAY  INTO  64  SUCCESSIVE  CONTROL  NENORY  WORDS. 

I 

ITO  KEEP  order  STRAIGHT  FIX  THE  BYTES  IN  THE  ARRAY  WITH 
ITHE  'ORDERS*  ROUTINE  BEFORE  STORING  WITH  THIS  ROUTINE, 

I 

FPI=  ARRAY  NUMBER 

FP2s  M'(8-BIT  field*  (MUST  BE  A MULTIPLE  OF 

OP=  address  of  first  CONTROL  MEMORY  WORD 

BAL.RO  STBYTS 


(ENTRY: 

( 

( 

( 

I 

(DISTURBS 

( 

t 


X,C*FL1«ASL*BL«FP2 
FP2  IS  INCREMENTED  BY 


(TIMING: 

134  MICROSEC  (526  NANOSEC/BYTE 

STORED) 

( 

ENTRY 

stbyts 

0000 

0000 

34AOOOOO 

stbyts 

LI 

BL,8 

60  THRU  OUTER  LOOP  TWICE 

0001 

OOCl 

30800009 

LRR 

ASL.DP 

SAVE  OP 

0002 

0002 

38014807 

GEN, 32 

X'38014807' 

PUT  AMRl  IN  BYTE  MOCE 

0003 

0003 

3F030014 

STBYTSl 

LOOP. 4 

STBYTS2 

0004 

0004 

47E08eA5 

L.W 

X,FP12« 

LOAD  XdNCR  FP2 

0005 

0005 

21C0A0FB 

LCU 

X(0) 

0006 

0006 

30020000 

SR 

C,0(OP) 

0007 

0007 

21C1A0FB 

LCU 

X(l) 

0000 

0008 

30020008 

SR 

C*e(DP) 

0009 

0009 

21C2A0FB 

LCU 

X(2) 

OOOA 

OOOA 

30020010 

SR 

C,16(DP) 

OOOB 

OOOB 

21C3A0FB 

LCU 

X<3» 

OOOC 

OOOC 

30020018 

SR 

C,24(0P) 

0000 

OOOD 

21C4A0FB 

LCU 

X(4> 

OOOE 

OOCE 

30020020 

SR 

C«32(0P) 

OOOF 

OOOF 

21C5A0FB 

LCU 

X(St 

OOlO 

OOlO 

30020028 

SR 

C*40(OP) 

0011 

0011 

21C6A0FB 

LCU 

X(6) 

0012 

0012 

30040030 

SR 

C«48(0P) <2 

INCR  DP 

0013 

0013 

21C7A0FB 

LCU 

X(7) 

0014 

0014 

30050037 

STBYTS2 

SR 

C,55(DP) t3 

INCR  DP.DECR  BL 

0015 

0015 

3040000A 

LRR 

OP.ASL 

RESTORE  OP 

0016 

0016 

29010019 

BZ.BL 

STBYTS3 

EXIT  IF  DONE 

0017 

0017 

28140001 

I NCR 

OP 

ELSE  INCR  OP 

ooia 

0018 

28010003 

B 

STBYTSl 

and  repeat  outer  loop 

0019 

0019 

380148FF 

STBYTS3 

GEN, 32 

X'380148FF* 

RESTORE  AMKl  TO  WORD 

OOlA 

OOIA 

28080000 

B 

0(R0) 

RETURN  TO  CALLER 

Figure  3-18.  Dump  Array  (STBYT$)  Listing 


The  line  thinning  algorithm  is  identical  to  that  described 
in  subsection  3.5  of  the  Final  Technical  Report  titled,  Associ a- 
tive  Array  Processing  of  Raster  Scanned  Data  for  Automated  Carto- 
graphy . ETL-0046.  However,  at  the  end  of  each  pass  through  the 
line  thinning  algorithm  (vertical  or  horizontal  array  sweep)  a 
pointer  CHNG  is  set  to  denote  at  least  one  cell  was  removed.  At 
the  end  of  four  passes  this  pointer  is  checked  and  the  THIN  routine 
is  exited  if  CHNG  is  zero. 

The  line  thinning  process  thins  in  such  a way  as  to  produce 
some  deformation  at  line  junctions  which  causes  indentations  in 
the  resulting  line  symbols  at  these  points.  This  is  noticeable 
at  4-mil  resolution,  but  is  much  less  apparent  at  2-mil  and  1-mil 
resolution.  The  reasoning  is  as  follows. 

A close  examination  of  the  digitized  source  data  showed  that 
some  line  junctions  produce  a 'build  up'  of  cells  after  the 
scanning/digitizing  process.  This  'build  up'  appears  to  be  related 
to  the  acuteness  of  the  angle  between  the  two  line  segments  and  is 
shown  below. 


A.  SOURCE  DATA  B.  SCANNED- REPLOTTED  SOURCE 
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IT' 


Because  the  line  thinning  algorithm  reduces  the  line  segments 
to  the  most  central  set  of  points  the  result  becomes  that  shown 
below. 


Generation  of  the  correct  line  symbology  from  this  data,  in 
this  case  a double  line,  produces  the  indented  effect  visible  at 
some  junctions  on  the  final  results  as  shown  in  the  following  diagram. 

Originally  techniques  were  investigated  to  remove  the 
indentations  by  either  breaking  the  lines  at  their  junctions  or 
massaging  the  junction  area. 

However,  these  investigations  didn't  produce  good  results. 

One  of  the  reasons  was  that  in  order  to  identify  a line  junction 
for  subsequent  segmenting  the  line  had  to  be  thinned  introducing 
the  indented  effect.  Secondly^  the  massaging  technique  involved 
thickening  and  thinning  which  again  produced  the  indentations. 

Because  it  is  not  possible  to  change  the  line  thinning  algorithm 
to  produce  a more  uniform  line  junction  it  appears  that  better 
results  are  obtained  by  using  one  or  more  of  the  following  techniques. 

0 Utilize  thinner  or  more  sharply  defined  source  material, 

0 scan/digitize  the  source  data  at  the  higher  resolutions,  or 

0 adjust  the  threshold  on  the  scanner/plotter  to  minimize 
the  addition  of  cells  at  line  junctions. 


i 

1 


,1 


j't 

: 1 
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3.9  VARIABLE  LINE  THICKENING  ROUTINES 

3.9.1  Functi on 

These  routines  generate  the  different  solid  line  symbols. 

3.9.2  Approach 

The  approach  was  to  generate  two  routines.  The  first 
routine  (THK)  performs  a series  of  4-cell,  2-cell,  and  1-cell  line 
thickness  similar  to  the  three  routines  described  in  subsection 
3.3.7  Line  Symbol  Generation  Routines  of  the  Final  Technical 
Report  ETL-CR-74-20  Associative  Array  Processing  for  Topographic 
Data  Reduction.  The  second  routine  (ELLSET)  generates  the  input 
parameters  to  THK.  The  input  parameters  provide  the  number  of 
passes  required  by  each  of  the  three  thickening  subroutines.  These 
parameters  are  determined  from  the  parameter  card  input  to  the 
Fortran  preprocessing  routine  PRPRC  which  indicates  the  thickness 
in  mils  of  the  final  line  symbol. 

The  three  line  symbol  patterns  used  to  perform  the  thickening 
are  shown  below. 


X 


1-CELL  2-CELL  4-CELL 

X 'ON'  CELL  OF  THINNED  LINE 

RESULTING  'ON'  CELLS  GENERATED 
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Figure  3-79  shows  the  top-level  flow  diagram  for  the  Variable 
line  thickening  Routine. 

3.10  DOUBLE  LINE  SYMBOL  ROUTINES 

3.10.1  Functi on 

These  routines  generate  the  double  line  symbols. 

3.10.2  Approach 

The  original  approach  to  Double  Line  Symbol  Generation  was 
to  generate  two  lines  within  the  arrays  ' on- the-f ly ' . The  larger 
width  line  corresponded  to  the  outside  diameter  of  the  required 
line  symbol,  and  the  smaller  line  was  the  width  of  the  inside  edges 
of  the  double  line  symbol.  For  a detailed  explanation  the  reader 
is  referenced  to  subsection  3. 3. 7. 3 of  the  Final  Technical  Report 
ETL-CR-74-20  Associative  Array  Processing  for  Topographic  Data 
Reducti on . 

Because  of  the  varying  degree  of  widths  (outer  and  inner 
diameters  of  the  double  line  symbols)  which  may  occur,  it  was 
decided  to  use  a different  technique.  The  technique  will  be  to 
make  two  separate  skid  line  symbols  whose  widths  correspond  to  the 
outer  and  inner  diameters  of  the  double  line  symbol.  This  will 
be  accomplished  by  using  the  Variable  Line  Thickening  routines 
discussed  in  subsection  3.9.  The  smaller  line  symbol  data  will 
then  be  exclusively  or'd  with  the  larger  line  symbol  data  at  merge 
time  and  thus  generate  the  required  double  line  symbol. 

3.11  LINE  EDITING  ROUTINES 
3.11.1  Functi on 

The  function  of  the  line  editing  routines  are  to  improve 
the  line  quality  of  data  prior  to  processing  by; 

1)  Joining  small  line  breaks,  and 

2)  Remove  'pips'  on  the  thinned  line  data. 
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3.11.2 


The  approach  was  twofold.  First,  the  line  joining  is 
accomplished  by  performing  a variable  line  thickening  on  thinned 
line  data,  and  then  rethinning.  This  function  accomplishes  the 
same  results  as  described  in  subsection  3.3.6  Line  Break  Detection 
and  Correction  of  the  Final  Technical  Report  ETL-CR-74-20  Associative 
Array  Processing  for  Topographic  Data  Reduction. 


Removal  of  line  'pips'  is  performed  by  detecting  line 
segments  less  than  some  predetermined  length  (based  on  resolution) 
and  removing  those  segments  by  vectorizing  the  line  segments  in 
question  and  deleting  the  points  along  the  line  as  they  are  found. 

3.11.3  Implementati on 

The  line  joining  routines  are  EDITSET  a level  zero  routine 
and  THKED  a level  3 routine.  These  routines  are  identical  to  the 
Variable  Line  Thickening  routines,  except  that  the  initial  para- 
meters for  THKED  are  obtained  from  a parameter  card  input  to  PRPRC 
(preprocessing  routine)  which  identifies  the  amount  of  edit 
thickening  required. 

The  'pip'  removal  routine  CLIP  is  a level  3 routine. 

3.12  BROKEN  LINE  SYMBOLOGY 

3.12.1  Functi on 

The  function  of  the  Broken  Line  Symbology  routines  is  to 
create  symbology  for  trails,  intermittent  streams  and  second 
class  roads.  Because  these  symbols  consist  of  a regular  pattern 
of  on/off  cells,  they  were  called  Broken  Line  Symbols. 

3.12.2  Approach 

The  approach  was  to  preserve  as  much  of  the  existing 
software  as  possible.  (Refer  to  ETL-0046  subsection  3.8  for 
details  of  existing  software.)  Because  creating  railroad  symbology 
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was,  in  reality,  somewhat  unrelated,  its  function  was  removed  and 
placed  in  a separate  module.  To  allow  for  more  dense  data,  the 
"setting  of  paint  bits"  section  of  the  module  was  modified  to 
operate  on  a block  basis.  This  modification  allows  array  vectors 
to  be  off-loaded  to  the  disk  at  the  end  of  a block.  In  the  old 
method,  all  array  vectors  were  required  to  be  core-resident  at  the 
time  paint  bits  were  set. 

3.12.3  Implement ati on 

Three  routines  are  required  to  create  Broken  Line  Symbols. 
They  are  CALRS,  STPTBT  and  PAINTl.  Each  module  is  discussed  in 
detail.  The  CALRS  calculate  the  master  vector  lengths  and  other 
running  sum  values  required  by  the  set  paint  bit  routine.  It  is 
called  only  once  for  each  sheet,  (after  vectorizing  and  prior  to 
STPTBT.)  It  performs  the  following  functions: 

1.  Retrieve  number  of  master  vectors  from  MASTAB+1 
and  store  in  NOMVEC. 

2.  Retrieve  number  of  array  vectors  from  MASTAB  and 
store  in  NOAVEC.  Load  array  vector  number  of 
increments  from  VNITAB  in  master  vector  order  until 
N0AVEC=0,  along  with  the  number  of  increments, 
save  the  array  vector  number  and  the  reverse  bit, 
RVBT,  in  the  array. 

3.  Flag  each  group  by  setting  the  first  bit  of  each 
group  in  START  flag. 

4.  Continue  steps  2 and  3 until  N0MVEC=0  or  there  is 
no  more  array  space. 

5.  Calculate  the  running  sum  values  for  each  Master 
vector  by  calling  the  routine  VLSUM$.  This  running 
sum  becomes  the  E,  or  end,  value  for  each  array 
vector  in  the  set  paint  bit  routine. 

6.  Replicate  the  last  value  in  each  group  (the  master 
vector  length  or  MVL)  by  calling  the  routine  REPVAL$. 

7.  Output  the  E,  MVL,  and  RVBT  to  VNITAB  using  the 
array  vector  number  to  address  VNITAB. 
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8.  If  there  are  more  master  vectors,  repeat  steps 
2-7,  else,  exit. 

NOTE 

If  a master  vector  contains  more  than  256  array 
vectors,  the  message  ERROR  AT  X'nnnn'  (where 
'nnnn'  is  an  address  in  CALRS)  will  be  printed 
out  on  the  terminal  and  the  program  will  halt. 

Figure  3-20  shows  the  vector  number  of  increments  table 
before  and  after  the  execution  of  CALRS. 

The  set  paint  bit  routine  (STPTBT)  is  used  to  create  the 
broken  line  symbology.  It  is  executed  on  a block  basis,  setting 
paint  bits  for  all  array  vectors  in  a particular  block.  The 
routine  CALRS  must  have  been  run  prior  to  the  execution  of  this 
routine.  The  array  vectors  (ARVECTAB)  for  a particular  block  and 
the  corresponding  address  table  (VADTAB)  must  have  been  read  off 
the  disk  prior  to  the  call. 

This  routine  is  similar  to  the  old  set  paint  bit  routine 
with  the  following  exceptions. 

1.  The  end  (E)  and  master  vector  length  (MVL)  values 
are  retrieved  from  VNITAB. 

2.  The  number  of  increments  (VI)  are  retrieved  from 
ARVECTAB  headers. 

3.  The  threshold  values  for  a particular  pattern  are 
variables  (rather  than  constants.)  They  are  stored 
in  a 5 word  table  which  is  inde.\ed  by  the  resolution 
of  the  data  being  created. 

The  PAINTl  routine  rasterizes  the  vectorized  data  back 


into  the  arrays  for  final  output.  It  is  identical  to  the  PAINT 
routine  developed  under  the  last  contract. 
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Figure  3-20.  VNITAB  Before  and  After  Execution  of  CALRS 


3.13  RAILROAD  SYMBOL  ROUTINES 


3.13.1  Functi on 

The  function  of  the  railroad  symbol  routines  is  to  create 
railroad  symbology  (tick  marks)  from  railroad  source  data.  The 
railroad  tick  marks  must  be  equally  spaced  and  lie  perpendicular 
to  the  railroad  center  line. 

3.13.2  Approach 

The  basic  approach  to  the  creation  of  railroad  symbols  is 
to  locate  the  position  of  each  tick  mark  along  the  center  line  and 
it's  angle  of  rotation.  This  data  is  then  used  as  input  to  the 
Point  Symbol  Routines  to  generate  the  tick  marks.  The  Railroad 
center  lines  will  be  generated  using  the  Variable  Line  Thickening 
Routines  discussed  in  subsection  3.9. 

Because  this  approach  involves  balancing  and  the  location 
of  specific  points  along  a line,  vectori zati on  is  employed  as  in 
broken  line  symbols. 

3.13.3  Impl ementati on 

3.13.3.1  General  - Creation  of  the  Railroad  tick  marks  involves 
the  following  steps: 

1.  Vectorize  input. 

2.  Link  to  from  Master  Vectors. 

3.  Balance  and  locate  the  tick  marks, 

4.  Compute  the  angle  of  rotation  of  the  tick  marks, 
and  their  X,Y  locations. 

5.  Draw  the  tick  marks  using  the  Point  Symbol  Routines. 
These  steps  are  discussed  in  more  detail  in  the  following  paragraphs. 

3.13.3.2  Vectori ze  - This  technique  is  similar  to  that  described 
in  subsection  3.16. 

3.13.3.3  Li nk  - This  routine  is  the  one  discussed  in  subsection 
3.16. 


3.13.3.4  Balance  - This  routine  is  a level  1 routine  called  TCKED 
(see  Figure  3-21)  which  finds  the  locations  of  the  railroad  ties 
(ticks)  in  a similar  manner  to  the  SETPTBT  routine  discussed  in 
subsections  3.12  and  3. 8. 3. 5 of  the  previous  Technical  Report, 
Associative  Array  Processing  of  Raster  Scanned  Data  for  Automated 
Cartography  (ETL-0046).  TCKFD  utilizes  the  results  of  the  CALRS 
routine  which  calculates  the  master  vector  lengths  and  the  running 
sum  values.  The  TCKFD  routine  is  similar  to  the  original  SETPTBT 
routine  (as  used  for  Railroad  Symbology),  except  for  the  following 
changes : 

1.  TCKFD  is  executed  on  a block  basis. 

2.  The  end  (E)  and  master  vecotr  lengths  (MVL)  are  retrieved 
from  VNITAB. 

3.  The  number  of  increments  in  each  array  vector  (VI)  are 
retrieved  from  the  ARVECTAB  headers. 

4.  The  railroad  symbol  threshold  values  are  variable  with 
respect  to  resolution,  as  are  the  ends  of  the  railroad 
ticks. 

5.  During  the  determination  of  the  ticks  the  end  values  LE, 
>RE  overwrite  the  array  vector  increments  within  the 
ARVECTAB,  and  the  number  of  ticks  and  their  start 
address  within  ARVECTAB  are  saved  in  VADTAB. 

6.  The  final  tick  values  (at  the  end  of  each  data  block) 
are  located  as  pairs  of  LE,  RE  values  stored  in  a tick 
location  table  TKLTAB  where  the  number  of  ticks  and  the 
initial  address  relative  to  the  start  of  TKLTAB  overwrite 
the  VNITAB. 

3.13.3.5  Compute  - This  routine  is  a Level  0 Routine  called 
(CRE8TK),  see  Figure  3-22,  which  generates  the  angle  of  rotation 
and  the  absolute  location  of  the  center  of  the  ticks. 

The  first  part  of  the  routine  is  similar  to  the  first  part 
of  the  Create  Tick  Routine  documented  in  Section  3.9.36  of  the 
previous  Technical  Report,  Associative  Array  Processing  of  Raster 
Scanned  Data  for  Automated  Cartography. 
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Figure  3-21.  Find  Tick  (TCKFD)  Routine  Flow  Diagram 
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Figure  3-21.  Find  Tick  (TCKFD)  Routine  Flow  Diagram 
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Figure  3-21.  Find  Tick  (TCKFD)  Routine  Flow  Diagram 
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Figure  3-22.  Create  tick  (CRE8TK)  Routine  Flow  Diagram 
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Figure  3-22.  Create  tick  (CRESTK)  Routine  Flow  Diagram 
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The  second  part  of  the  routine  generates  the  angle  of 
rotation  of  the  ticks  and  the  absolute  location  of  their  centers. 

The  absolute  value  of  the  center  of  each  tick  is  determined 

2 2 '^D 

where  (X^.Y^)  and  2)  are  the  absolute  locations  of  the  two 

end  points  of  each  tick  from  the  top  left  hand  corner  of  the  sheet. 
The  table  lookup  values  to  approximate  the  angle  of  rotation  of 
the  ticks,  are  determined  from  the  tables  shown  in  Figure  3-23. 

This  provides  an  accuracy  of  ^7.5  degrees.  The  angle  of  rotation 
is  computed  from  Xjj/Yp=0  and  the  true  angle  of  rotation  is; 

90+9  for  Xjj/Yq  negative  valued 
90-0  for  Xq/Yjj  positive  valued 

The  angle  of  rotation  and  the  absolute  locations  of  the  centers 
of  the  railroad  ticks  are  output  to  a buffer  in  Bulk  Core  Memory 
which  corresponds  to  the  buffer  used  in  Prepro  (see  subsection  3.15). 
A message  is  then  sent  to  the  user  to  load  and  execute  the  Point 
Symbol  Routines. 

3. 14  AREA  FILL  ROUTINES 

3.14.1  Functi on 

The  area  fill  routines  are  routines  which  fill  in  outlined 
areas  defined  by  'ON'  cells  with  a solid  mass  of  'ON'  cells. 

3.14.2  Approach 

The  approach  was  to  utilize  the  existing  area  fill  routines 
described  in  subsection  3.10  of  the  Final  Technical  Report  ETL-0046 
Associative  Array  Processing  of  Raster  Scanned  Data  for  Automated 
Cartography.  However,  the  routines  which  generate  the  area  fill 
patterns  and  logically  .AND.  the  patterns  into  the  filled  areas 
have  been  excluded.  The  reason  for  this  is  as  follows. 


AO-A049  69a 


UNCLASSIFIED 


AO 

A049ese 


6000YEAR  AEROSPACE  CORP  AKRON  OhIO  F/O  9/5 

ASSOCIATIVE  ARRAY  PROCESSINO  OF  RASTER  SCANNED  DATA  FOR  AUTOMAT— ETC (U) 
NOV  77  N J ADAMS*  J M VOCARr  K LOSCH  DAAa53-76^-0ia6 


aER-I6S23 


ETL-0132 


A.  Relationship  between  TanO/0 

B.  Range  for  each  0 


Figure  3-23.  Table  Lookup  Definition 

The  old  approach  of  generating  the  area  fill  patterns  from  some 
library  was  considered  inappropriate  for  a semi -producti on  system 
for  the  following  reasons. 

1.  Each  different  pattern  had  to  be  designed  as  a binary 
coded  image  and  saved  in  a Library.  The  design  of 
each  symbol  was  extremely  time  consuming,  particularly 
at  higher  resolutions  and  yet  must  closely  adhere  to 
the  specifications  required  to  satisfy  the  Carto- 
graphic community. 

2.  If  the  pattern  size  exceeds  the  array  size  then 
subsections  of  the  pattern  had  to  be  generated  to 
match  across  the  array  boundaries. 

3.  In  the  case  of  patterns  which  consist  of  randomly 
located  sub-patterns  such  as  swamps,  marshland,  etc. 
this  library  storage  technique,  and  the  calling  of 
each  sub-pattern  into  the  arrays  when  required, 
became  extremely  complicated. 
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In  the  new  approach,  a map  sheet  sized  STIC-PAT  of  each  area  fill 
pattern  will  be  scanned  at  the  proper  resolution  and  a raster 
plot  tape  generated.  This  tape  will  then  be  merged  with  the 
output  tape  from  the  STARAN  as  an  off-line  procedure.  Where 
scanning  a full  sheet  requires  more  than  one  tape,  a smaller 
segment  of  the  pattern  will  be  scanned  and  then  merged  a multiple 
number  of  times  to  produce  the  correct  area  fill  pattern. 

This  technique  will  provide  the  following  advantages. 

1.  Each  area  fill  STIC-PAT  pattern  will  only  have  to  be 
scanned  once  for  each  required  resolution  and  then 
the  file  can  be  saved  and  used  when  needed. 

2.  The  quality  of  the  symbols  produced  will  be  more 
acceptable  to  the  Cartographic  community. 

3.  In  the  case  of  complicated  patterns,  the  randomness  of 
sub-patterns  (such  as  in  swamps)  will  be  as  good  as 
the  design  of  the  original  STIC-PAT. 


4.  New  patterns  can  be  added  to  the  area  fill  Library 
with  a minimal  of  design  time  and  subsequent  software 
effort. 

5.  As  a production  environment  is  approached  these  area 
fill  tapes  may  be  stored  as  resident  files  on  some 
larger  and  faster  peripheral  device  and  merged  with  the 
STARAN  output  tape  as  it  is  generated  on  a record- 
by-record  basis. 

3.14.3  Impl ementati on 

The  area  fill  routines  are  all  level  3 routines  and  are; 

1.  AFIL, 

2.  GTL,  and 

3.  SVL. 

These  routines  are  the  same  as  the  area  fill,  get  last  slice, 


-i 


.1 


and  save  last  slice  routines  described 
the  Final  Technical  Report  ETL-0046. 


in  subsection  3.10  of 
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3.15  POINT  SYMBOLOGY  ROUTINES 

3.15.1  Function 


The  function  of  these  routines  is  to  provide  a method  for 
creating  point  symbols  such  as,  churches,  houses,  schools,  route 
markers,  mine  symbols,  numbers,  etc.  at  the  proper  locations  on 
the  map.  The  churches,  schools,  houses  must  be  rotatable  and  all 
symbols  must  reflect  the  quality  and  accuracy  of  set  symbol 
standards  at  the  resolution  being  processed. 

3.15.2  Approach 

The  approach  was  to  modify  and  expand  the  existing  point 
symbology  routines  in  such  a manner  as  to  minimize  the  size  of  the 
point  symbols  defined  in  a Point  Symbology  Library  and  yet  obtain 
the  desired  results  at  the  resolution  required. 

3.15.3  Implementation 

3.15.3.1  General  - The  locations  of  the  point  symbols  to  be 
produced  on  a map  overlay  sheet  are  defined  by  an  X,Y  pair  and  an 
angle  of  rotation.  These  X,Y  values  are  in  inches  from  the  upper 
left  hand  corner  of  the  overlay;  the  angle  is  specified  in  integer 
degrees  between  0 and  359.  This  data  is  input  on  cards.  The 
actual  point  symbols  are  stored  as  aX,aY  pairs  from  an  arbitrarily 
defined  center  point.  These  points  are  contained  in  a Point 
Symbol  Library  and  describe  only  selected  points  of  each  symbol. 

This  technique  utilizes  the  same  set  of  library  points  to 
generate  the  different  resolutions  of  the  same  symbol  and  thus 
minimize,  disk  space  utilized,  access  time,  and  Library  generation 
complexi ty. 
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The  1-mil  symbols  are  created  by  rotating  and  translating 
these  points  to  correspond  to  their  correct  position  in  an  array 
and  data  block.  These  data  points  contained  in  the  array  are 
then  thickened  and  smoothed  to  achieve  a 1-mil  symbol  as  shown  in 
Figure  3-24. 

The  2-mil  symbol  is  created  from  the  same  set  of  library 
points,  but  the  values  of  those  points  are  divided  by  2 (a  shift 
of  1 bit). 

Thickening  is  utilized  to  achieve  the  2-mil  symbology  as 
shown  in  Figure  3-25. 

The  4-mil  symbols  are  created  from  the  same  set  of  points, 
as  shown  in  Figure  3-26,  but  the  output  points  are  divided  by  4 
(shift  of  2 bits). 

Some  of  the  points  used  in  generating  the  4-mil  symbols  are 
redundant  and  map  into  neighboring  locations  after  the  division 
by  4. 

Each  symbol  is  defined  with  a minimum  number  of  points  while 
retaining  the  appearance  of  the  correct  resolution. 

Figure  3-27  shows  the  Point  Symbol  editing/thickening  flow 
diagram.  The  thickening  subroutines  (Thick  1,  Thick  3)  are 
identical  to  the  Variable  line  thickening  routines  discussed  in 
subsection  3.9. 

The  1-mil  edit  subroutine  smoothes  out  the  boundaries  of  the 
Point  Symbols  to  create  the  1 mil  resolution  boundary.  Essentially, 
this  subroutine  generates  cells  (points)  when  the  boundary  changes 
direction  by  90°  and  both  of  the  perpendicular  segments  are  greater 
than  two  cells  long.  Figure  3-28  shows  an  example  of  this  editing 
process. 
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Figure  3-25,  2-Mil 
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The  aX,aY  points  are  rotated  the  amount  specified  on  the  input 
card  and  translated  to  the  X,Y  center.  The  rotated,  translated 
points  are  then  drawn  into  the  arrays  at  the  proper  time,  (based 
on  the  order  of  the  raster  output). 

3.15.3.2  Read  Input  Cards  - This  is  a Fortran  routine  INPTA  which 
reads  the  input  card  data.  The  form  of  the  input  cards  is  shown 
in  Figure  3-29.  This  routine  reads  the  X and  Y offset  for  all 

the  data.  (This  offset  is  used  to  align  point  symbols  on  scanned 
road  data).  It  then  reads  the  X and  Y coordinate  of  the  center  of 
the  symbols,  the  angle  of  rotation  (0  for  non-rotated  symbols),  and 
the  symbol  flags.  The  X and  Y coordinates  are  offset  and  converted 
to  STARAN  format.  The  converted  data  is  written  on  the  disk 
(F0R001.DAT)  for  use  by  the  Point  Symbol  Preprocessing  routine 
(PREPRO). 

3.15.3.3  Preprocess  and  Sort  Symbol  cells  - PREPRO  is  a STARAN- 
APPLE  Level  0 program.  It  is  used  to  read  the  disk  file  created 
by  INPTA  and  to  pre-process  all  the  point  symbol  inputs.  It  calls 
the  following  subroutines:  RDISKl  and  SINCOS. 

3.15.3.3.1  PREPRO  - This  routine  preprocesses  the  input  data  for 
the  point  symbol  routine.  It  calls  RDISKl  to  read  the  disk  file 
created  by  INPTA.  The  actual  (x,y)  values  of  the  centers  are 
converted  to  working  coordinates  (Block  and  column.  Array  and  row) 
pairs.  The  data  are  then  sorted  by  block  plus  overlap. 

The  basic  process  performed  in  this  routine  is  not  unlike  the 
process  performed  in  the  last  contract.  For  a detailed  explanation 
of  point  symbol  processing  refer  to  ETL-0046  Section  3.11.  The 
major  discepancies  are  outlined  below: 

1.  The  routine  now  handles  l-,2-  and  4-mil  resolutions. 
This  affects  the  conversion  from  x,y  to  working 
coordinates,  as  the  number  of  resolution  units  per 
inch  is  now  a variable. 
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2.  The  start  column  and  start  row  of  an  array  load, 
along  with  the  number  of  processed  columns  and  rows 
have  been  changed  to  variables  as  opposed  to 
constants.  These  values  are  used  in  the  conversion 
process.  They  were  changed  to  variables  to  allow 
for  changes  in  overlap. 

3.  The  sine  and  cosine  values  are  now  calculated  in  the 
STA?>AN  (as  opposed  to  the  PDP)  using  the  routine 
SINCOS. 

4.  In  order  to  allow  for  more  symbol  types,  the  symbol 
flag  field,  SYMBPT,  was  expanded  to  12  bits. 

5.  In  order  to  accommodate  more  paint  symbols  in  the 
future,  the  remaining  process  (rotation/translation) 
was  split  out  to  execute  on  a block  basis.  PREPRO 
must  sort  the  input  data  into  block  groups  to 
facilitate  this  change.  If  the  center  of  a symbol 
lies  in  an  overlap  region,  for  example,  between 
Block  2 and  Block  3,  the  information  pertaining  to 
that  symbol  will  be  duplicated  in  the  block  group 
for  2 and  for  3. 

6.  PREPRO  creates  BLKTAB  which  is  an  8 word  table  used 
by  BLKPRO.  A bit  set  in  that  table  indicates  that 
the  corresponding  block  has  data  to  be  processed. 

For  example.  Bit  3 in  word  0 of  BLKTAB  set  implies 
that  Block  3 contains  data  and  Bit  2 in  word  1 of 
BLKTAB  set  implies  that  Block  34  contains  data. 

3.15.3.3.2  SINCOS  - Is  a subroutine  called  by  PREPRO  which 
computes  the  sine,  cosine  of  the  angles  associated  with  the  symbols 
read  by  PREPRO. 

3.15.3.3.3  RDISKl  - Is  a subroutine  called  by  PREPRO  which  reads 
the  disk  file  F0R001.DAT  created  by  INPTA  into  a Bulk  Core  Memory 
buffer  whose  initial  address  is  defined  by  the  name  BUFFER.  It 
reads  124  bytes  at  a time;  120  of  which  contain  input  point  data. 

It  reads  until  an  EOF  is  found.  The  program  then  returns  control 
to  PREPRO.  (see  Figure  3-3o)- 
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3.15.3.4  Block  Processing  Point  Symbol  Routine  - BLKPRO  Is  a 
STARAN-APPLE  Level  1 program.  This  routine  processes  point  symbol 
data  on  a block  basis.  It  assumes  that  the  center  Information  has 
already  been  converted,  sorted,  and  stored  In  Bulk  Core  by  the 
routine  PREPRO.  The  main  functions  performed  by  this  routine  are: 

1.  Load  arrays  with  center  Information, 

2.  Search  on  symbol  type, 

3.  Bring  symbol  library  off  disk, 

4.  Rotate  Symbol  Points,  and 

5.  Sort  on  array  load  basis. 

3.15.3.4.1  BLKPRO  - BLKPRO  first  checks  the  appropriate  bit  in 
BLKTHB  (discussed  in  the  previous  subsection)  to  determine  if  there 
is  data  to  be  processed  this  block. 

If  there  is  data  to  be  processed,  the  BLOCK,  COL,  ARRAY,  ROW, 
SYMBPT,  SIN  and  COS  for  each  symbol  ^re  loaded  Into  the  arrays, 
one  symbol  per  array  word. 

A search  is  then  made  on  the  symbol  flag,  SYMBPT.  When  a 
symbol  is  found,  RDISKl  is  called  to  read  the  appropriate  symbol 
file  from  the  disk. 

GETSYM  is  then  called  to  load  the  symbol  data  into  the  arrays. 

PROSYM  is  called  to  rotate  and  translate  the  symbol  points 
and  to  store  them  into  bulk  core. 

If  a symbol  is  larger  than  128  points,  the  calls  to  GETSYM 
andPROSYM  are  repeated  until  the  entire  symbol  has  been  processed. 

The  search  on  SYMBPT  is  repeated  until  all  symbol  types  have 
either  been  checked  and/or  processed. 

The  output  data  is  then  ordered  for  drawing.  The  ordering 
technique  is  the  same  as  described  under  the  previous  effort. 
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3.15.3.4.2  PROSY’  - PROSYM  rotates  (If  the  symbol  is  rotable)  and 
translates  the  symbol  points.  It  also  outputs  the  processed  data 
to  a buffer  in  bulk  core. 

PROSYM  varies  from  the  old  routine  in  two  respects; 

1.  The  data  can  be  4-, 2-  or  1-mil.  Therefore,  prior  to 
output  the  data  is  adjusted  to  the  appropriate 
resolution.  The  4-mil  data  is  divided  by  4 (a  2 bit 
shift)  the  2-mil  data  is  divided  by  2 (1  bit  shift) 
and  the  1 mil  data  is  unchanged.  The  divide  is 
performed  at  the  end  of  the  processing,  so  that 

the  fractional  portion  of  the  aX,aY  can  participate 
in  the  rotation  and  translation  process. 

2.  Because  the  creation  of  2 and  1 mil  symbols  ulti- 
mately requires  thickening,  data  must  be  drawn  in 
the  overlap  region  of  the  array  to  insure  proper 
symbolization.  Therefore,  the  sort  routine  was 
modified  to  include  points  in  the  overlap  region. 

When  PROSYM  is  exited,  all  symbols  of  the  type  being  processed 
have  been  rotated  and  translated.  (That  Is,  if  there  are  more 
symbol  centers  of  the  type  or  than  can  be  processed  in  1 pass, 

(1  pass=4  symbols  at  GAC,  8 at  ETL)  PROSYM  will  loop  until  all  are 
done. ) 

3.15.3.4.3  GETSYM  - The  GETSYM  routine  is  identical  to  the 
GETSYM  routine  of  the  last  effort  with  one  exception.  Symbol  size 
has  been  increased  to  128  points,  (as  opposed  to  64)  so  that  each 
array  can  hold  only  2 sets  of  symbol  points.  GETSYM  loads  the 
symbol  point  data  from  bulk  core  into  the  first  128  words  of  each 
available  array  and  then  replicates  the  data  by  rotating  and  .OR.ing. 
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3.15.3.5  Draw  Symbol  Routine  - The  DRAW  routine  is  identical  to 
that  routine  created  under  the  last  contract  effort.  It  was  the 
array  table  ARYTAB  created  by  BLKPRO  to  determine  whether  there 

is  any  data  to  be  drawn  in  the  array.  (ARYTAB  acts  just  like  BLKT 
BLKTAB  except  it  indicates  an  array  load  of  data  instead  of  a 
Block.)  Data  is  drawn  by  setting  a bit  in  the  array  for  each 
COL,  ROW  pair  stored  in  the  bulk  core  buffer  by  PROSYM. 

3.15.3.6  Edit  Point  Symbols  - PTEDT  - This  routine  is  a Level  3 
routine  called  PTEDT  which  generates  the  final  point  symbols  from 
the  set  of  aX.AY  values  positioned  within  the  arrays  by  DRAW  and 
edits  this  data  in  different  ways  depending  on  the  resolution. 

This  editing  is  required  to  fill  in  holes  resulting  in  some  aX.aY 
cells  mapping  into  each  other  during  the  rotation  and  drawing 
stages.  The  top-level  flow  diagram  is  shown  in  Figure  3-27. 

Each  of  the  three  point  symbol  generation  editing  subroutines 
performs  a slightly  different  function  depending  upon  the  resolu- 
ti  on . 

3. Is. 3. 6.1  Four-Mil  Edit  - This  subroutine  performs  the  same 
editing  process  that  was  described  in  subsection  3.14.3.5  of  the 
last  Technical  Report. 

3.15.3.6.2  Two-Mil  Edit  - This  subroutine  performs  the  following 
edit  process. 

1.  Perform  a 3 x 3 cell  thickening  process  on  the  symbols 
within  the  arrays. 

2.  Complement  the  data  within  the  arrays  (ON  cells  become 
OFF  cells,  OFF  cells  become  ON  cells). 

3.  Perform  a 3 x 3 cell  thickening  process  on  the  comple- 
mented symbols. 

4.  Complement  the  data  within  the  arrays. 

Figure  3-31  shows  the  results  of  the  4 steps. 
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Figure  3-31.  Symbol  Editing  process  for  2-mi1,  1-mil  data. 


3. lb. 3. 6. 3 One-Mil  Edit  - This  subroutine  performs  the  same 
editing  process  as  the  ;^-mil  edit  except  that  a 5 x 5 cell  thick- 
ening pattern  is  used  instead  of  the  3x3  cell  pattern.  A 
second  editing  function  is  required  to  smooth  the  edges  of  the 
one-mil  symbols.  Essentially,  this  subroutine  improves  the 
resolution  of  the  symbols.  Figure  3-32  shows  the  flow  diagram 
of  this  subroutine.  The  routine  locates  changes  in  direction  along 
the  symbols  boundaries.  These  directional  changes  must  be  by 
90®  and  continue  (in  both  directions)  for  at  least  two  cells.  The 
point  nearest  to  the  corner  is  then  set,  thus  resulting  in  a 
smoothed  corner.  Steps  1-3  are  the  Line  Separation  Routine.  It 
is  this  routine  which  has  been  rewritten  in  the  following  manner. 

The  two-cell  cross  reduce  CRSRDC,  and  one-cell  box  reduce 
BXRDC  routines  are  used  together  with  their  corresponding  two  cell 
and  one  cell  restore  routines  CRSRST,  BXRST.  However,  the 
'reduced'  data  is  completely  unloaded  from  the  arrays  and  saved  in 
an  area  of  Bulk  Core  Memory.  After  further  reductions,  to  eliminate 
the  line  weight  desired,  this  data  is  then  exclusive  OR'd  back 
into  the  arrays  to  produce  the  desired  line  weight. 

Besides  the  parameter  to  isolate  the  desired  line  weight 
WGTDES,  a second  parameter  is  input.  This  parameter  FATFLG 
indicates  that  the  line  weight  desired  corresponds  to  the  thickest 
line  on  the  sheet.  Using  this  second  parameter  reduces  the  steps 
involved  in  separation  when  isolating  the  largest  line.  Figure  3-33 
shows  the  Line  Separation  Flow  Diagram. 

3.16  LINE  SEPARATION  BY  THICKNESS  ROUTINES 

3.16.1  Functi on 

The  function  of  the  Line  Separation  routines  is  to  take  one 
raster  input  tape  and  extract  a desired  line  weight  selected  prior 
to  processing  by  the  user.  The  extracted  line  is  represented  by 
its  enter  line  and  stored  in  raster  from  on  an  output  tape. 
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32.  Improve  Resolution  of  1-mil  Symbols  Flow  Diagram 
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-33.  Line  Separation  Routine  Flow  Diagram 
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3.16.2  Approach 

The  approach  was  to  modify  the  existing  Line  Separation 
Process  described  in  Section  3.12  of  the  Technical  Report  entitled 
Associative  Array  Processing  of  Raster  Scanned  Data  for  Automated 
Cartography. 

3.16.3  Impl ementati on 

The  line  separation  technique  involves. 

1.  Reducing  all  line  widths  until  lines  smaller  than 
the  line  weight  desired  are  eliminated. 

2.  Reducing  the  remaining  lines  until  the  desired 
line  is  eliminated. 

3.  "Subtracting"  the  lines  remaining  in  Step  2 from 
the  set  produced  in  Step  1.  The  result  is  the  line 
weight  desired. 

4.  Editing  the  line  by  joining  line  breaks  and 
eliminating  clutter. 

3.17  VECTORIZATION 

3.17.1  Functi on 

The  function  of  vectorizing  is  to  convert  raster  data  to 
vector  data. 

3.17.2  Approach 

The  approach  to  this  problem  was  to  utilize  the  routines 
written  during  the  last  contract  DAAK02-75-C-0114  with  modifications 
to  allow  the  processing  of  4-mil,  2-mil,  and  1-mil  data. 

3.17.3  Implementation 

This  subsection  describes  the  changes  to  the  routines 
documented  in  subsection  3.13  of  the  previous  Technical  Report, 
Associative  Array  Processing  of  Raster  Scanned  Data  for  Automated 
Cartography.  The  routines  previously  written  were: 
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1.  Boundary  and  End  Point  Tagging  Routine. 

2.  Line  Following  Routine. 

3.  Vector  Link  Table  Update  Routine. 

4.  Master  Vector  List  Creation  Routine. 

3.17.3.1  Array  Boundary  and  End  Point  Tagging  Routine  - This 
routine  has  now  been  split  Into  two  Level  3 routines  to  facilitate 
the  execution  of  these  routines  utilizing  the  alternating-paging 
approach  discussed  In  subsection  3.3  (Executive  Routine)  of  this 
report.  The  resulting  two  routines  are: 

a)  Junction/Detection  and  End  Point  Tagging  Routine  (OTCTA) 

b)  Array  Boundary  Tagging  Routine  (DTCB) 

3.17.3.1.1  DTCTA  - This  routine  processes  the  thinned  line  data. 

The  routine  detects  and  removes  a 3 x 3 set  of  cells  at  all 
junction  points  and  then  tags  and  identifies,  by  position,  all 
end  points.  A detailed  flow  diagram  is  shown  in  Figure  3-34; 
also,  a detailed  description  follows. 

An  initial  check  is  made  to  determine  if  the  line  data  in 
the  array(s)  coincides  with  data  along  any  edge  of  the  overlay 
being  processed.  In  this  case  slices  of  data  (equivalent  to  half 
of  the  overlap  between  data  blocks)  occurring  along  the  correspon- 
ding edge  are  cleared.  This  is  to  ensure  that  any  array  vectors 
that  cross  these  boundaries  are  not  tagged  as  array  boundary  vectors. 

The  next  step  is  to  clear  a frame  around  the  data  to  be 
processed,  namely  word  slices  0,1,2,3,252,253,254,255,  and  bit 
slices  0,1,  — 31.  This  provides  working  storage  and  locations 
for  the  results. 

Two  bit  slices  (254,255)  are  now  set  up  to  denote  those  array 
words  which  will  participate  in  the  junction  point  or  end  point 
detection  process. 


Figure  3-34.  Detailed  Flow  Diagram  of  DICTA 

(sheet  1 of  3) 
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Figure  3-34. 
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The  criteria  for  the  existence  of  a junction  point  is  that: 

1.  The  center  point  is  'ON'  (i.e.,  all). 

2.  The  sum  of  the  surrounding  eight  points  is  three  or  four. 

3.  The  sum  of  the  'logical  OR'  between  each  pair  of 
surrounding  points  is  six. 


The  next  step  is  to  locate  and  tag  all  the  end  points  by 
their  word  and  bit  slice  positions.  Again,  the  process  is  the 
same  as  for  'line  editing'  except  that  the  number  of  surrounding 
points  is  one  (the  sum  of  the  logical  OR's  of  each  pair  is  auto- 
matically two),  as  shown  in  Figure  3-35.  While  the  end  points 
are  being  located  any  'spots'  which  have  no  surrounding  points 
are  edited  out  of  the  array. 


WORD  POSITION  OF 
END  POINT  WITHIN 
BIT  SLICE 


Figure  3-35.  End  Point  of  Line 


The  bit  slice  and  word  positions  of  each  end  point  are  then 
stored  in  consecutive  array  words  in  the  field  positions  shown  in 
Figure  3-36.  Each  end  point  is  also  tagged  by  setting  its  corres 
^ ponding  bit  in  the  end  point  flag  slice  (EM)  and  end  point  or 

array  boundary  vector  slice  (MM).  These  flag  slices  are  used  by 
the  line  following  routine.  The  resulting  array  layout  of  the 
data  after  tagging  is  shown  in  Figure  3-37. 
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Figure  3-36.  Array  Layout  After  Tagging 


3.17.3.1.2  DTCTB  - This  routine  processes  the  data  resulting  from 
the  junction  point  detection  and  end  point  tagging  routine  (DTCTA). 

The  routine  locates  and  tags  all  the  array  boundary  vectors. 

The  steps  involved  in  locating  array  boundary  vectors  are 
given  below: 

1.  Take  an  exterior  boundary  slice  (top,  bottom,  left,  or 
right)  at  the  edge  of  the  data  to  be  vectorized  by  the 
Line  Following  Routine. 

2.  Locate  any  connected  array  vectors  on  the  interior 
inside  of  this  boundary  slice. 

3.  Store  the  locations  (FPl,  FP2)  of  any  such  connections 
below  the  end  point  values  determined  by  DTCTA,  and  in 
the  same  field. 

4.  Set  the  corresponding  tag  slice  MM  (slice  252),  and 
the  combined  end  point  and  array  boundary  tagged, 
slice  indicator. 

Steps  1 through  4 are  repeated  for  the  remaining  three  boundary 
slices.  Figure  3-38  shows  an  example  of  this  array  boundary 
testing  and  tagging. 

3.17.3.2  Line  Following  Routine  - This  routine  is  a Level  3 
routine  called  (LINF)  which  converts  the  raster  data  to  vector 
data  by  following  and  vectorizing  the  individual  line  segments 
with  the  arrays. 

The  routine  is  identical  to  the  routine  described  in  subsection 
3.13.3  of  the  Final  Technical  Report,  Associative  Array  Processing 
of  Raster  Scanned  Data  for  Automated  Cartography.  However,  the 
location  of  the  end  point  and  array  boundary  point  flag  slice  (MM) 
has  been  changed  from  slice  36  to  slice  252. 
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3.17.3.3  General  Table  Update  Routine  - This  routine  is  a Level  3 
routine  called  (TAB)  which  updates  the  Vertical  and  Horizontal 
Array  Boundary  Tables  on  an  iteration  basis.  The  Array  Vector 
Link  Table  is  then  updated  at  the  end  of  each  data  block  from 
these  boundary  tables.  Figure  3-39  shows  the  top-level  flow 
diagram  of  this  routine. 

The  base  addresses  of  the  vertical  and  horizontal  array  bound- 
aries are  initialized  during  the  first  i terati on/bl ock . There  are 
two  vertical  array  boundary  tables  and  one  horizontal  array 
boundary  table.  As  each  block  is  processed  the  vertical  array 
boundary  table  that  is  used  to  update  the  VLNKTAB  is  alternated 
between  the  two  available  tables.  These  three  tables  are  located 
in  Bulk  Core  Memory  as  shown  in  Figure  3-40.  All  array  vectors 
which  cross  art*y  boundaries  are  stored  into  the  next  available 
locations  in  the  array  boundary  tables.  After  the  last  two  arrays 
in  the  block  have  been  processed  the  routine  updates  the  Array 
Vector  Link  Table. 

Following  is  a detailed  description  of  the  Update  Routine,  as 
shown  in  Figure  3-41.  A mask  slice  MSKERST  is  initialized  to  reset 
bit  zero  and  bit  sixteen  of  those  array  vectors  which  are  linked 
to  another  array  vector.  By  this  method  an  end  of  a string  of 
connected  array  vectors  may  be  determined  (see  Master  Vector 
Creation  Routine).  The  locations  denoted  by  VL,  VR,  HL,  HR  are 
used  to  point  to  the  next  available  16-bit  locations  in  the 
Vertical  and  Horizontal  Array  Boundary  Tables  (where  L denotes  the 
left  side  and  R the  right  side).  These  are  initialized  at  the 
base  addresses  of  the  appropriate  array  boundary  tables. 

Figure  3.42  shows  an  example  of  the  use  of  these  tables  after 
the  first  and  second  blocks  have  been  processed.  Notice  that  VL 
and  VR  alternate  between  the  two  vertical  array  boundary  tables. 
After  the  first  block  has  been  processed  only  the  horizontal  array 
boundary  table  is  used  to  update  the  Link  table.  The  array  boundary 
tables  are  updated  from  the  array  vector  numbers  which  were  created 
during  the  Line  Following  Routine. 
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Figure  3-41.  Flow  Diagram  of  General  Table  Update 
Routine  (sheet  1 of  3) 
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Figure  3-41.  Flow  Diagram  of  General  Table  Update 
Routine  (sheet  2 of  3) 
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Figure  3-41. 


Flov#  Diagram  of  General  Table  Update 
Routine  (sheet  3 of  3) 
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Figure  3-43  shows  an  array  prior  to  this  step.  The  array 
boundary  table  together  with  the  sixteen  bit  half-word  that  will 
contain  the  array  vector  number  is  determined  by  the  four  array 
boundary  flag  slices  denoted  by  L (left),  T (top),  B (bottom), 
and  R (right). 

If  this  is  the  last  iteration  in  the  data  block  the  array 
vector  link  table  (VLNKTAB)  is  updated  from  the  boundary  tables. 

The  T/F  flag  bits  denoted  by  bit  positions  0 and  16  determine  how 
the  VLNKTAB  is  updated.  Figure  3-44  shows  an  example  of  the  four 
possible  combinations  of  T/F  values.  The  arrows  on  the  array 
vectors  A and  B denote  the  direction  in  which  they  were  vectorized. 

Figure  3-45  shows  how  the  VLNKTAB  is  updated  for  the  four 
different  cases.  Prior  to  storing  the  array  vectors  into  VLNKTAB 
the  T/F  flag  bits  are  'switched  off  by  utilizing  the  mask 
generated  at  the  beginning  of  the  routine.  The  steps  are  shown 
below: 

1.  Store  array  boundary  table  word  in  X response  store. 

2.  Generate  the  mask  in  Y response  store. 

3.  'Logically  AND'  X with  Y. 

4.  Output  X bits  0 through  31. 

Figure  3-46  shows  this  process.  When  all  words  within  the  array 
boundary  tables  have  been  used  to  update  the  VLNKTAB  the  routine 
is  completed. 

3.17.3.4  Master  Vector  Table  Routine  - This  routine  is  a Level  0 
routine  called  MVECT  which  generates  the  master  vector  list  from 
the  Array  Vector  Link  Table.  Figure  3-47  shows  the  top-level  flow 
diagram  of  the  Master  Vector  Table  Routine.  Following  is  a detailed 
description  of  the  new  routine.  Figure  3-48  shows  the  associated 
Flow  Diagram. 
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Figure  3-47,  Top  Level  Diagram  of  Master  Vector  List  Routine 

128 


W 


NXTVC 


A.VEC  AN 
^END  PT^ 


DOES  THE^ 
A. VECTOR 
REQUIRE 
^EVERSINa- 
\ 


ANOTHER  A.VEC 
v LINKED  TO  . 
\J^IS  ONE/^ 


NWMHDR 

SET  NEW  HEADER 


IT^AN 
END  PT 


CLIND 


Figure  3-48.  Flow  Diagram  for  Master  Vector  Routine  (sheet  2 of  8) 
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Flow  Diagram  for  Master  Vector  Routine  (sheet  4 of  8) 


Figure  3-48.  Flow  Diagram  for  Master  Vector  Routine  (sheet  5 of  8) 
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Figure  3-48.  Flow  Diagram  for  Master  Vector  Routine  (sheet  7 of  8) 
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The  operations  involved  in  the  generation  of  the  Master  Vector 
Table  are  as  follows. 

1.  Initialize  the  number  of  array  vectors  generated  (NPT). 
Then,  set  all  bits  in  array  0 to  one.  These  bits  are 
used  on  a slice-by-slice  basis  to  identify  which  array 
vectors  have  not  been  output  as  part  of  a master  vector. 

Z.  Locate  the  next  array  vector  not  output  already  to 
master  vector  list. 

3.  Follow  the  array  vector  linkage  reversing  array  vectors 
where  necessary  (see  Figure  3-49)  until  either: 

a)  An  end  point  is  located 

b)  The  starting  array  vector  in  a closed  loop  string 
is  hit  for  the  second  time  (see  Figure  3-50). 

4.  Output  all  the  array  vectors  in  this  string  to  the  master 
vector  table,  and  identified  as  a unique  master  vector 
number. 

5.  Delete  all  the  associated  array  vector  flags  from 
array  0. 

6.  Repeat  steps  2 through  5 until  all  array  vectors  have 
been  output  (NPT  = zero). 

7.  Update  the  header  of  the  master  vector  table  with  the 
resolution,  and  the  number  of  master  vectors  generated. 

3.18  CONTOUR  TAGGING  INVESTIGATION 
3. 18. 2 Functi on 

The  function  of  this  investigation  was  to  determine  tech- 
niques for  tagging  vectorized  contours  generated  from  raster 
scanned  contour  maps. 

3.18.2  General 


Figure  3-51  depicts  the  general  data  flow  conceived  for 
contour  to  grid  process  (CONTAGRID).  The  concepts  and  techniques 
that  have  been  developed  to  perform  this  program  are  fully  docu- 
mented in  the  PROPOSAL  FOR  CONTOUR  TAGGING  AND  GRIDDING  USING 
STARAN/CDC  SOFTWARE.  Many  of  the  conclusions  reached  are  a direct 
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Contour  overlays  will  be  scanned  on  either  the  DMA  or  ETL 
scanner  plotters  at  2-mil  resolution.  The  resulting  tapes  will  be 
read  in  on  a block  basis  by  the  CDC  and  broken  up  into  four  array- 
load groups.  These  data  will  then  be  stored  on  the  844  disk. 

As  needed,  the  array  load  groups  will  be  read  from  the  disk 
and  transferred  via  the  command  channel  interface  to  the  STARAN. 

Here,  the  data  will  be  loaded  into  the  arrays  (converted  from  run- 
length  to  binary)  and  the  contour  tagging  process  will  begin. 

The  first  part  of  the  Contagrid  process  requires  identification 
of  index  contours  as  opposed  to  non-index  contours.  This  will  be 
done  by  separation  based  on  line  weight. 

After  separation  the  data  will  be  vectorized.  Array  vector 
tables  and  address  tables  will  be  shipped  to  CDC  core  and  stored 
on  the  844  disk  as  buffer  space  becomes  filled.  Vector  headers 
will  reflect  the  index  or  non-index  status  of  the  line.  Data 
required  to  form  master  vectors  and  the  master  vectors  themselves 
will  be  stored  on  the  STARAN  RK05  disk. 

Once  the  master  vectors  have  been  formed,  the  automatic  edit- 
ing processes  will  begin.  Numerics  will  be  identified  and  flagged 
for  deletion.  Depression  will  be  recognized,  ticks  flagged  for 
deletion  and  intermediate  vectors  joined  to  form  continuous  contours. 
Deep  cuts  will  be  recognized  and  joined  appropriately  to  corres- 
ponding master  vectors.  Breaks  caused  by  the  placement  of  numerics 
on  the  contours  will  be  located.  Appropriate  endpoints  will  be 
paired  and  joined  in  a straight  line  manner.  The  joining  array 
vectors  will  be  shipped  to  the  CDC  for  disk  storage  and  new  master 
vectors  will  be  created  which  reflect  this  linkage.  Should  regis- 
tration marks  along  the  sheet  boundaries  be  required  by  the  gridding 
process,  they  will  be  identified  and  flagged  at  this  time. 

The  contour  data  will  now  be  ready  for  manual  editing.  Data 
will  be  formatted  for  output  to  plotting  devices.  CalComp  format- 
ting will  be  performed  by  the  CDC.  (Any  raster  plotting  required 
will  be  formatted  by  the  STARAN.)  Numerics,  flagging  areas  of  change 
(adds  or  deletes)  and  problem  areas  (open  endpoints),  will  be  used 
to  communicate  corrections  back  into  the  system. 
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Next,  the  elevation  data  will  be  read  in  and  associated  with 
the  open  index  contours.  As  part  of  the  manual  input,  a list  will 
be  made  (counter  clockwise  from  the  upper  left-hand  corner  of  the 
sheet)  of  the  elevations  of  all  index  contours  which  cross  the 
sheet  boundaries.  The  open  contours  will  be  ordered  in  a counter- 
clockwise manner.  By  matching  the  input  list  (card  format)  to  the 
master  vectors  flagged  as  indices  in  the  ordered  list  in  the  STARAN, 
the  association  of  elevation  to  linear  list  can  be  easily  made. 

When  the  open  contours  (opens)  have  been  appropriately  tagged, 
the  closed  contours  will  be  tagged  by  finding  the  order  of  embedded 
contours,  and  their  proximity  to  known  opens.  This  process  will  involve 
rasterization  of  contour  data.  Prior  to  rasterization,  the  master 
vector  numbers  must  be  associated  with  the  array  vectors  which  com- 
prise the  master  vectors.  Because  the  array  vectors  are  resident  on 
the  CDC  disk,  the  association  will  be  made  in  the  CDC.  The  tagging 
of  closed  contours  will  then  proceed. 

If,  after  the  contours  have  been  automatically  tagged,  there 
remain  contours  with  unknown  elevations,  these  will  be  formatted 
for  output  to  the  CalComp  where  the  editor  will  assign  appropriate 
elevations.  The  data  will  then  be  prepared  for  the  gridding  process. 
The  tagged  contour  vector  strings  will  be  reduced  to  the  resolution 
of  the  grid  output  (10-mils).  Saddle  and  ridge  (S/R)  lines  and 
spot  elevations  will  be  input  via  CDC  tape  and  converted  from  DGR 
format  to  a format  compatible  with  the  contour  data.  Intersections 
of  S/R  lines  with  contour  lines  will  be  found  and  a grid  of  known 
elevation  points  will  be  created.  (It  should  be  noted  that  the 
modularity  of  the  Contagrid  software  will  allow  for  direct  input  of 
DGR  tagged  contour  data  into  the  gridding  portion  of  the  system, 
bypassing  the  auto-tagging  software  entirely.)  This  grid  of  points 
will  be  partitioned  into  array  load  groups  and  shipped  to  the 
STARAN  where  planar  interpolation  will  be  used  to  find  the  missing 
elevation  points. 

It  is  highly  unlikely  that  the  input  sheet  will  be  so  poorly 
registered  on  the  scanner  as  to  require  alignment  correction. 
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However,  since  skew  could  be  a problem  area  a possible  solution  is 
addressed  at  the  end  of  the  software  subsystems  section.  If  such 
corrections  are  necessary,  they  will  be  performed  prior  to  the 
tagging . 

The  final  grid  will  be  sent  back  to  the  CDC  and  formatted  into 
DGR  format.  The  final  output  will  be  written  onto  magnetic  tape. 

3.18.3  Line  Separation  of  index  contours  - It  is  intended  to  use 
the  Line  Separation  routines  discussed  in  subsection  3,16  to  perform 
the  separation  of  the  index  and  non-index  contours. 

In  order  to  determine  the  reaction  of  these  routines  to  2-mil 
contour  data  the  following  has  been  performed. 

Goodyear  Aerospace  has  separated  a 2-mil  data  tape  of  contours 
generated  on  the  DMA  scanner/plotter.  The  results  of  this  process 
are  being  forwarded  to  ETL  for  plotting  on  the  scanner/plotter. 
Conclusions  will  be  drawn  from  the  resulting  plots. 
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APPENDIX  A - MAGNETIC  TAPE  FORMATS 

1.  GENERAL 

This  appendix  describes  the  tape  record  formats  used  to 
provide  the  input/output  data  to  the  STARAN  AAP-raster  pro- 
cessing. 

There  are  four  Input/Output  forms  that  can  be  accommodated 
by  the  current  software;  these  are: 

1)  ETL  Plot  Format  Input  data  - ETL  Plot  Format  - 
Output  data, 

2)  ETL  Run-Length  Coded  (RLC)  Input  data  - ETL  Plot 
format  Output  data, 

3)  DMA  Run-Length  Coded  Input  data  - DMA  Run-Length 
Coded  Output  data,  and 

4)  DMA  Run-Length  Coded  Input  data  - ETL  Plot  Output 
data. 

Each  of  the  three  formats  - 1)  ETL  Run-Length  Coded  Format,  2)  ETL 
Plotter  Format,  and  3)  DMA  Run-Length  Coded  Format  - are  described 
below. 

2.  ETL  SCANNER  RLC  FORMAT 

Figure  A-1  shows  the  ETL  RLC  format.  The  input  data  consists 
of  four  8-bit  bytes  per  transition  from  bl ack- to-whi te  or  white- 
to-black.  Each  byte  corresponds  to  one  frame  of  a nine  track 
magnetic  tape.  The  X dimension  denotes  a circumferential  spot 
position  on  any  track  and  the  Y dimension  corresponds  to  the 
track  (scan-line).  A 'positional  marker'  is  generated  every 
250  spot  positions.  This  is  identical  to  a transition  point 
and  the  color  code  corresponds  to  the  color  of  the  previous  transi- 
tion. 

3.  ETL  PLOTTER  FORMAT 

Plot  tape  format  uses  1 byte  (8-bits  + parity)  of  data  to 
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control  8 light  emitting  diodes  for  simultaneous  exposure  of  8 
tracks  on  the  drum.  Thus,  one  tape  record  contains  the  data  for 
all  circumferential  positions  of  8 successive  tracks.  As  with 
binary  scanning,  the  plot  records  control  each  circumferential 
spot.  The  bits  controlling  the  light  emitting  diodes  are  either 
set  (1)  or  not  (0),  indicating  exposure  or  not. 

Records  lengths  are  computed  as  with  binary  scanning,  by 
the  number  of  spot  positions  around  the  drum  (for  one  track). 

The  number  of  records  is  equal  to  1/8  the  number  of  scan  records, 
since  8 tracks  are  exposed  by  each  record. 
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Figure  A-2.  ETL  Plotter  Format 

4.  DMA  SCANNER/PLOTTER  RLC  FORMAT 

Figure  A-3  shows  the  DMA  RLC  format.  The  Run-Length  value 
(1  - 127^g)  represents  the  number  of  spot  increments  to  power-on 
or  power-off  the  laser  beam  along  the  X-coordinate  (scan)  line. 
Each  Magnetic  Tape  (MT)  RLC  record  contains  a 16-bit  Y-address, 
a mode  bit,  the  word  count,  and  the  RLC  data.  The  Y-address 
represents  a scan  line.  The  mode  is  0 for  RLC  data.  The  word 
count  defines  the  actual  number  of  16-bit  data  words  which  follow 
the  word  count. 
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Figure  A-3.  DMA  Run-Length-Coded  (RLC)  Formet 
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Two  magnetic  tapes  are  generated  by  the  DMA  scanner/plotter. 
One  tape  contains  data  with  even  Y-addresses,  and  one  contains 
data  with  odd  Y-addresses. 
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APPENDIX  B - I/O  MERGING/SEPARATING  ROUTINES 


The  STARAN  S-500  system  supports  only  two  tape  drives; 
therefore,  two  off-line  utility  routines  were  written  to  handle 
DMA  run-length  coded  data  tapes.  The  functions  of  these  two 
routines  were: 

1)  to  merge  two  DMA  run-length-coded  input  tapes 
into  one  input  tape  containing  alternate  records 
from  each  tape,  and 

2)  to  separate  a STARAN  processed  DMA  run-length-coded 

output  tape  into  two  tapes,  one  containing  records 
(0,2,4,...)  etc.,  and  the  other  containing  records 
(1,3,5 ) etc. 

These  utility  routines  are  written  in  FORTRAN  and  are  executed 
on  the  XEROX  SIGMA  9 which  currently  supports  three  tape  drives. 


APPENDIX  C - COMTAL  DEBUG  AID 


The  Comtal  Debug  Aid  routine  (COMTAL)  is  an  APPLE  program 
executed  from  Bulk  Core  Memory.  This  routine  is  called  by  the 
Executive  when  requested.  It  transfers  a pair  of  arrays  of  data 
to  the  left-half  of  the  Comtal  screen,  and  subsequently  an 
adjacent  pair  of  arrays  of  data  to  the  right-half  of  the  screen. 
This  data  is  in  raster-binary  format  and  is  generated  as  a 
graphic  overlay  (black  (off),  color  (on))  on  the  screen. 

The  user  generates  requests  for  viewable  data  via  a list  of 
four  parameters.  These  parameters  identify: 

1)  Pass  Number  (1,2,)  etc., 

2)  Data  Block  Number  (0,1,2,...)  etc., 

3)  Iteration  Number  within  data  block  (0,1,2,...)  and 

4)  Function  within  iteration  (1,2,3...)  etc.,  to  be 
displayed . 

The  parameters  are  input  via  the  keyboard  (SDM)  as  32-bit  values 
to  the  HSDB,  where  each  8-bits  of  each  word  identifies  one  of  the 
four  parameters.  The  first  HSDB  location  is  identified  by  the 
label  RCHK  and  is  at  location  X'680'^g. 

Exampl e 

Display  the  results  after  line  thinning,  and  after 
line  symbol  generation  during  the  processing  of  track  data. 
Display  these  results  for  the  first  and  second  iterations 
of  the  first  data  block.  The  resulting  parameters  generated 
in  successive  words  of  HSDB  would  be: 


Location 
X' 680' 

X'68r 
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pass 


data  block 


0 0 
0 1 


0 0 
0 0 


0 0 
0 0 


0 2 
0 2 


i terati on 
function 
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The  current  parameters  pass  (p),  data  block  (d),  iteration 
(i),  and  function(f)  are  constantly  monitored  by  the  executive. 
When  the  above  parameters  match  those  currently  being  checked  in 
the  HSDB  list,  then  the  data  in  the  arrays  is  output  to  the 
Comtal  screen  (left-half).  Later,  during  processing,  the  adjacent 
two  arrays  of  data  (p,  d+1,  i,f)  are  output  to  the  screen  (right- 
half),  to  generate  a 2 x 2 set  of  arrays  on  the  screen. 

The  Comtal  contains  three  graphic  displays,  green,  blue,  and 
red.  The  viewable  data  sets  (four  arrays)  are  output  to  the 
graphics  alternating  from  green  to  blue,  and  then  to  red.  As 
each  colored  graphic  is  displayed,  the  preceding  graphic  is 
saved  and  may  be  retrieved  later  by  the  user.  However,  once  all 
three  graphics  have  been  used,  the  next  request  for  graphic 
data  will  overwrite  the  first  of  the  three  graphics  saved,  etc. 
Thus,  only  three  graphic  data  sets  can  be  saved  for  viewing  at 
any  one  time.  The  parameter  list  must  also  be  generated  in 
ascending  value  with  respect  to  the  numbers  associated  with  p, 
d,  i,  and  f. 

This  routine  allows  the  user  to: 


1)  view  the  results  of  some  Raster  Processing  function. 
Thinning,  Symbol,  Generation,  etc.,  on-line  without 
processing  the  whole  overlay  sheet  and  obtaining 
some  hard  copy  output, 

2)  terminate  a processing  procedure  if  a problem  is 
encountered  due  to: 

a)  bug  in  a program, 

b)  strange  data  form  encountered,  and 

3)  demonstrate  various  aspects  of  Raster  Processing. 
While  this  routine  is  not  directly  usable  by  ETL  personnel,  it 
has  been  saved  as  an  example  of  a possible  technique  for  future 
implementation  in  the  ETL  Facility. 
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